X-581-73-369 


9** 



70 $35 


(t* 

|Q0 

1 un 


& 

! :s 


r 


i 

l 


A 


w 

(0 00 
r— I lO 

U m 

rt sfr 

O CN 


AE-C ATTITUDE DETERMINATION AND 
CONTROL PRELAUNCH ANALYSIS 
AND OPERATIONS PLAN 


rn 

X 

rn 

o 



R. D. WERKING 
R. D. HEADRICK 
C. F. MANDERS 
R. D. WOOLLEY 



GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 



RECEDING PAGE 


BLAN£ not filmed 


FOREWORD 


This document has been prepared by the Attitude Determination and Control 
Section and Computer Sciences Corporation to provide a description of attitude 
control support being supplied by the Mission and Data Operations Directorate. 
Included are descriptions of the computer programs being used to support the 
missions for attitude determination, prediction, and control. In addition, de- 
scriptions of the operating procedures which will be used to accomplish mis- 
sion objectives are provided. 
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SECTION 1 - INTRODUCTION 


1. 1 MISSION PLAN 

The objective of Atmosphere Explorer-C (AE-C) is to continue investigation of 
the upper atmosphere begun by AE-1 and AE-2, earlier spacecraft. It is de- 
signed to determine energy input at lower altitudes (120 to 300 kilometers) by 
measuring structural parameters (neutral and charged particle composition, 
density, and temperature), solar ultraviolet input, and airglow line intensities. 

The AE-C orbit has an initial perigee of 157 kilometers, an apogee of 4300 kil- 
ometers, and an inclination of 68. 1 degrees, chosen to give a slow rotation 
(-0.61 degree per day) of the argument of perigee. The latest available launch 
elements are listed in Appendix A. Perigee will be lowered periodically to 
altitudes as low as 120 kilometers for 1-day periods. These low-altitude oper- 
ations are possible because of the onboard Orbit Adjust Propulsion System 
(OAPS) and the additional stability provided by the momentum wheel. To assure 
performance of the mission, the spacecraft attitude must be maintained within 
the following constraints: 

• The spin-axis attitude is to be maintained within 2 degrees (arc) of 
the orbit normal (either positive or negative to maintain the Z-axis 
Sun angle less than or equal to 90 degrees). 

• The total angular momentum must be maintained within 10 percent 
of the nominal value of 1200 in. -lb-sec. 

• These attitude constraints will be met by operation of MSAD/AE 
for attitude determination, primarily on near real-time (apogee) 
passes and by Multi- Satellite Attitude Prediction/Atmosphere 
Explorer (MSAP/AE) for attitude control command generations. 
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1.2 ATTITUDE SENSORS AND CONTROL SYSTEMS 


Attitude sensing is provided by wheel horizon scanners (WHSs); a body horizon 
scanner (BHS), as shown in Figure 1-1; a fan-shaped digital solar aspect sen- 
sor; a magnetometer triad; and an elevation gimbal. Data from the Miniature 
Electrostatic Accelerometer (MESA) will be provided to the Orbit Determination 
group during OAPS firings. 

The horizon sensors respond to infrared radiation in the 14- to-16-micron CO^ 
band, chosen to minimize cloud noise. The WHS detectors receive radiation 
which is reflected off mirrors rotating on the momentum wheel at effective 
polar mounting angles of 104 and 118 degrees with respect to the spin (Z) axis. 
The WHS telemetry consists of Earth widths (elapsed time) for both sensors, 
split-to -index time (equivalent to a rotational angle), and wheel spin period. 

The BHS sensor has a similar detector mounted at a 110-degree polar angle to 
the Z-axis. A second BHS sensor is mounted at the same angle for redundancy, 
but data from only one of them is telemetered down in the form of Earth width 
and spin period. 

The Digital Solar Aspect Indicator (DSAI) has a fan-shaped field of view (FOV) 
which covers Sun angles, /3 , with 1-degree resolution from 0 to 180 degrees 
from he spin axis. The solar gate sensor enables the DSAI when the Sun is 
within 0.75 degree arc of the DSAI Y-Z plane. 

A set of biaxial gimbals will be provided to aim a solar experiment to the center 
of the Sun with an accuracy of 1 minute of arc. The elevation gimbal angle 
with respect to the spacecraft will be determined with an accuracy of ±0.5 de- 
gree. 

An orthogonal triad of magnetometers is planned for the AE-C mission: one 
will be mounted parallel to the +Z-axis and the others will lie in the X-Y plane, 
90 degrees out of phase. Each component will be sampled in telemetry at the 
main frame rate. No specification is made against magnetic materials, and 
spacecraft magnetic fields may be substantial. 
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An orthogonal triad of MESAs will be available to measure acceleration due to 

thrusting and atmospheric drag. One will be mounted along the Z-axis and the 

other two will have their sensitive axes at ±45 degrees to the X-axis. There 

will be three constraintment ranges (A, B, and C), the most sensitive of which 

-5 8 

(C) permits measurements from 2 x 10 g (full scale) to 4x10 g, Con- 

_3 

straintment range A permits measurement from 8 x 10 g (full scale) to 
-5 

1.6 x 10 g . This range will normally be used during spacecraft thruster 
firing. 

The AE spacecraft has two modes of operation: the pitch loop control mode 
and a freewheeling, or tachometer, loop mode. In the pitch loop mode, the 
spacecraft rotates about the spin axis (Z), at either a 4— revolution— per— minute 
or despun (Earth-pointing) mode of 1 revolution per orbit. In the tachometer 
mode, the spin rate can be selected by ground command from 1 rpm to 
approximately 8 rpm. In the despun mode, the spin axis is normal to the 
orbital plane and the spacecraft Y -axis can be aligned with the local vertical 
or at any commanded pitch angle from this position. 

The attitude control system consists of a momentum wheel assembly, magnetic 
torquers, nutation dampers, automatic roll control, pitch control, and yaw 
thruster. 

Figure 1-1 shows the location of the attitude control system components. 

Momentum consists of the sum of the body and the wheel inertia-velocity prod- 
ucts. The simplest and most accurate (±0. 1 percent) determination of total 
momentum can be made from the telemetered wheel period when despun. Dur- 
ing contact passes while spinning, both wheel and body rate data are also avail- 
able from the horizon and Sun sensors. Thus, a running record of body inertia 
can be kept, so that subsequent momentum determination during the spin mode 
is possible. 
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The momentum wheel provides a momentum reference for stabilizing the 
spacecraft and a reaction torque for controlling the spin rate and despun atti- 
tude. The initial angular momentum will be attained by spinning the spacecraft, 
using the delta vehicle. Momentum will be maintained by appropriate transfers 
between the spacecraft and the wheel. (In the despun mode, essentially all the 
momentum will be transferred to the wheel, but in the spin mode, it will be 
distributed between the wheel and the spacecraft. ) The spin rate will be set to 
the rate commanded from the ground by sensing the rotation rate with respect 
to the local vertical and controlling wheel speed appropriately. 

The magnetic torquers control spacecraft orientation and momentum magnitude. 
After horizon sensors have determined the orientation, a magnetic precession 
torque, generated by a Z-axis coil either automatically or by ground command, 
is applied to the spin axis to align the axis normal to the orbit plane. Momentum 
magnitude control will be accomplished via a magnetic torque applied about the 
spin axis. Appropriate switching of the magnetic field of a spin coil in the 
spacecraft X-Y plane provides an interaction with the Earth's magnetic field 
to effect the desired momentum magnitude change. 

Passive nutation dampers on the despun portion of the spacecraft will be used 
to damp spacecraft oscillations induced by aerodynamic drag, thruster firing, 
and other disturbing forces. 

During low-perigee passes and apogee restoration maneuvers, the primary at- 
titude disturbance results in roll deviations in the vicinity of perigee. Automatic 
roll control is provided for backup in case of missed station contacts to reduce 
this attitude error by automatic sensing and magnetic torquing. Because roll 
deviations occurring at perigee appear as roll deviations at the subsequent 
apogee, closed-loop (i. e. , automatic) roll sensing will take place at apogee. 

Spin or pitch angle control is provided in either a closed-loop mode or in an 
open-loop (or tachometer) mode. In the closed-loop mode, control is achieved 
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by momentum transfer between the flywheel and the main body by means of a 
direct current torque motor, in the open-loop mode, the desired value of the 
wheel speed can be set by ground command at any of 32 values from 10 to 41 
radians per second, to allow control of the body spin rate in the range of 1 to 
8 rpm. The closed-loop mode is employed either during the despun mode (one 
revolution per orbit) and at a spin rate of 4 rpm. 

Pitch error is measured by one of the two wheel horizon sensors. The sky- 
Earth and Earth-sky transitions are used to generate a pulse that bisects the 
two horizon pulses. This pulse represents a local vertical reference that can 
then be compared with the encoder-generated index pulse to obtain the pitch 
error. 

The necessary horizon splitting is accomplished by gating a clock to a counter 
for sensor Earth-time duration and, beginning with the next sky- Earth transi- 
tion, counting down at twice the frequency. This process produces a time- 
bisecting pulse during the next wheel spin. 

Two redundant pitch axis control loops will be available. Either of the torque 
motors with its corresponding electronic loop can be selected by ground com- 
mand and either horizon sensor can be used with the selected loop. 

The pitch electromechanical assembly for pitch control consists of an outer 
fixed housing that interfaces with the spacecraft and of an inner rotating turn- 
table. Mounted onto these two items are an encoder with redundant coils, 
labyrinth seals, redundant brushless torque motors; an oil reservoir; the fly- 
wheel; a mirror assembly; and bolometers and associated electronics. The 
turntable assembly is supported by the external housing through a single bearing. 

1.3 GROUND DATA FLOW 

The ground data flow is shown in Figure 1-2. 
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Figure 1-2. Ground Data Flow 


ORBIT-CHANGE AND YAW- MANEUVER REQUESTS 





SECTION 2 - MISSION ANALYSIS 


To maintain the attitude and momentum within mission constraints, the estab- 
lished procedure is followed: attitude determination, prediction, and control. 
Because of the relatively strong torques and tight constraints, special care 
must be taken in attitude accuracy, bias determination, and estimation of space- 
craft parameters. 

2.1 ATTITUDE DETERMINATION 

The basic method used for real-time spin-axis attitude determination is the 
conic intersection technique for both deterministic and least -squares process- 
ing. Sensor data from WHS, Sun sensors, and magnetometers are converted 
into cone angles about the reference vectors. Special filtering methods are 
used to average the cone angles to reduce the sensitivity to body dynamics. In 
terms of vectors on a unit sphere, the equations can be written as: 

Z • E = cos n 

Z . S = cos j3 (2-1) 

/*v 

Z • B = cos X 

/\ 

where Z = spacecraft spin-axis unit vector 
/\ 

E = satellite-Earth unit vector 

A 

S = satellite-Sun unit vector 
B = Earth's B-field unit vector 
n = nadir cone angle between Z and E 
p = Sun angle between Z and S 

/N /\ 

X = magnetometer cone angle between Z and B 

The horizon scanners, which operate in the 14-to-16-micron band, respond to 
radiation from the CO molecules in the upper atmosphere, and thus are not 
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subject to terminator problems that are common to optical sensors. They are 
sensitive, however, to other error sources, such as threshold drift, radiance 
variations, orbit error, and Earth oblateness. The Earth widths from the WHS 
telemetry are converted to nadir angles from the spherical triangle relation- 
ship: 


cos p = cos y cos n + sin y sin n cos f (2-2) 

c 2 2 

where y = sensor mounting (cant) angle, 118 degrees 

f = half width of the Earth from f = 77 EW2/SP 
Z 2 

EW2 = measured Earth width from WHS2 
SP = selected spin period 

P = angular radius of the Earth, calculated from the following: 


R + HT 


CO. 


sin p = 


where r = distance from the center of the Earth to the satellite 

HTqq = height of the CO effective horizon 
2 * 

Rg = radius of the Earth at the averaged horizon-crossing points 
The correction has been made for oblateness using a two-parameter model: 



e i 3tn *LAT 


+ C 2 Sln *LAT 


where R E q = equatorial radius 

= flattening coefficient 
0 lat = horizon -crossing latitude 

e 2 = empirical coefficient for seasonal effects 
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Although Equation (2-2) yields ambiguous values for n , the two solutions are 
easily resolved when the attitude is near the orbit normal. With both WHS sen- 
sors operational, the term cos p can be eliminated, removing the dependence 
on orbit ephemeris. The resulting dual-scanner equation is 

cos 7 9 - cos 7 1 

tan n = — — — ■=■ : — — — 

sin y cos f - sin y cos f 

JL Li imi 

without ambiguity. At apogee, the orbit and oblateness errors are small and 
the WHS2 single- scanner solution is more accurate (o^ =0.3 degree) because 
of fewer sensor-electronic errors. At apogee, the WHS2 sensor can tolerate 
roll angles of -9 degrees (n = 81 degrees) before missing the Earth. 

Sun angles are read directly from the DSAI Gray-coded telemetry with 1-degree 
resolution or an accuracy of 0^ = 0. 5 degree. The DSAI range is from 0 to 
180 degrees, and thus will provide data as long as the spacecraft is in sunlight 
and in the spinning mode. When the SPS is tracking the Sun, the elevation 
gimbal gives a separate Sun-angle measurement to about the same accuracy, 
but one which is more subject to performance errors. The gimbal range is 
0 to 90 degrees and can provide data also when the body is despun. 

The magnetometer cone angle, X , is obtained from 


cos X = 


H 

z 


/ 


2 2 2 

H + H + H 
x y z 


These magnetometers are subject to rather large errors due to quantization 

alone (cr = 1.25 degrees), and spacecraft magnetic fields are expected to be 
X 

large and variable due to coils, transmitters, tape recorders, experiments, 
and the iron momentum wheel. Constant biases can be resolved and removed 
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by routines in the telemetry processor (BAMN and RESIDU). The magnetom- 
eters will be used as backup data tor another reference vector for conic inter- 
section when the spacecraft is in shadow. Although the magnetometer data is 
available at the main-frame rate (16 per second), it will ordinarily be sifted 
down to about 1 per second to reduce the data volume. 

2.2 DISTURBANCE TORQUES 

A nominal launch will place the spin axis at negative orbit normal, and to 
satisfy mission requirements, the spin axis must track the orbit normal to 
within 2 degrees. The orbit normal describes a cone in inertial space, with the 
right ascension decreasing about 1. 5 degrees per day because of the Earth's 
oblateness. 


In addition to the tracking error caused by attitude sensing and control limita- 
tions, there are several external torques which precess the spin axis away from 
its proper orientation. 

The torques considered include aerodynamic, magnetic, solar pressure, grav- 
ity gradient, and the delta-V thruster. The effects of these disturbance torques 
are assessed by use of MSAP/AE simulation runs. 


The simulations were made with an orbit generator using the following prelim- 
inary orbit elements : 


Epoch 

Epoch time 

Semimajor axis 

Eccentricity 

Inclination 

Mean anomaly 

Argument of perigee 

Right ascension of ascending node 


December 7, 1973 
5. 53 hours U.T. 

8625. 1195 kilometers 

0.2421825 

68.1104 

359.4129 degrees 
165. 7542 degrees 
236.8955 


and the initial attitude at negative orbit normal. 
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One-day simulation runs were made 3 days after launch date by turning on 
various torque-calculation indicators in the MSAP satellite-dependent input. 
The effect of the varying perigee height was studied by changing the semimajor 
axis and the eccentricity to achieve the desired perigee altitude. The results 
are presented by MSAP/AE-C attitude summary and various MSAP plots. The 
summary of the attitude drift rate due to these torques is given in Table 2-1. 
An example of atypical run (120 kilometers, aerodynamic torque) is given in 
Appendix B. The increase in momentum in this run was caused by having the 
integration stepsize a multiple of the spin period. 

Because of the low perigee altitude, aerodynamic torque will be the primary 
cause of attitude disturbance for the AE-C mission. 

The aerodynamic torque is modeled by the following equation, using the simple 
drag-force model: 


T 

aero 


= r x F 


and 



/\ 


(- V) 


where 

P 

V 

A 

P 

/N 

V 
r 


drag coefficient =2.0 

atmospheric density 

magnitude of the satellite velocity 

surface area of the satellite projected on the plane perpendicular 
to V 

unit vector in the direction of satellite velocity 

the vector displacement of the satellite center of pressure from 
the center of mass 
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Table 2-1. Attitude Drift Rate Summary 


PERIGEE 

HEIGHT 

(KM) 

TYPE OF TORQUE 

SPIN RATE 
{RPM) 

CONDITIONS ! 

4a 

(DEG) 

A6 

(DEG) 

AL 

(IN.-LB -SEC) 

120 

AERODYNAMIC 

4 

C Q =2. YCMOFF =0.1" 

4.9 

1,4 

57.7 



0 

ZCMOFF = 0.1" 

5.0 

1.4 

98.9 

161 

AERODYNAMIC 

4 

C D =2. YCMOFF =0.1" 

0.4 

0.1 

1.7 



0 

ZCMOFF =0.1" 

0.4 

0.1 

8.6 

161 

MAGNETIC 

4 

ZB IAS = 500 POLE-CM 

0.1 

0.0 

0 



4 

ZBIAS = 1000 POLE-CM 

0.2 

0.0 

0 | 

161 

MAGNETIC 

0 

XBIAS =50 POLE-CM 

0.0 

0.0 

0 

161 

GRAVITY 

4 

NEGATIVE ORBIT NORMAL 

0.0 

0.0 

0 




2-DEGREE CONING 

0.0 

0.0 

0.1 

161 

SOLAR 

4 

ZCMOFF =0.1" 

0.0 

0.0 

0 

161 

NONE 

4 


0.0 

0.0 

0 




It is assumed here that shadowing effects which might be caused by solar point- 
ing subsystem are negligible. The projected area of the AE-C satellite surface 
is given by 


7TD 

A = I hD sin 0. | +-— — cos 6. 
p i 4 i 

where D = diameter of cylinder 
h = length of cylinder 
0. = angle between V and cylinder axis 

The aerodynamic force is proportional to the atmospheric density at perigee. 
The drift rate is proportional to the force and the offset of the centroid from the 
center of mass. 

The other disturbance torque which is observable from the drift rates is the 
magnetic torque, given by 


T = M x B 
mag 

where M = magnetic dipole moment 
B = Earth's magnetic field 

Cases were simulated for residual dipole moments of 500 and 1000 pole-cm 
along the Z-axis, with results as shown in Table 2-1. In the despun mode, a 
dipole moment of 50 pole-cm along the X-axis produced no effect. 

When the delta-V thruster is fired it will produce a disturbance torque 

T =~r x? 


2-7 



where F = force vector along the thrust axis 

r = thrust offset vector from the center of mass of the spacecraft 

Because F is nearly along the X-axis, the center of mass offset along x will 
not produce any torque. The torque from the y and z components can be 
evaluated separately for roll and spin-rate changes. 

Using the impulse approximation, 

AH = Az • F At 
y x 

When the firing occurs along or against the velocity vector of 

AH 

A( P r = ±h L 

TOTAL 

at perigee or apogee, this results in a change in roll angle, <p . 

For typical values (Az = 0.1 inch, F^ = 4.0 pounds, At = 260 seconds, and 
^TOTAL = in. -lb-sec) , we have 5.0 degrees . For perigee- 

lowering maneuvers, it will be necessary to preroll the spacecraft to avoid 
entering the atmosphere with the attitude far from the orbit normal. 

The Ay offset gives a spin rate change 

AH = - Ay F At 
z 

of -104.0 in. -lb-sec for Ay = 0.1 inch. 

2.3 CONTROL TORQUES 

An optimal magnetic attitude control method has been developed for AE-C to 
maintain the spin-axis direction within 2 degrees of orbit normal. The control 
law has been shown to be optimal in the sense of minimum time, minimum 
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energy, and minimum attitude deviation. The deviation and some eomparitive 
results is presented in the paper included as Appendix C. 

The yaw thruster is provided for emergency maneuvers, and produces a torque 

T =~r XF 

where F - force vector along the thrust axis 

~r - thrust offset vector from the center of mass of the spacecraft 

The thruster is located on the spacecraft to produce a strong torque along the 
Y-axis that can be used to perform a 180-degree yaw maneuver and can also be 
used for rapid correction of the spin-axis orientation. 

The Magnetic Attitude Spin Coil (MASC) is commutated to produce a torque along 
the Z-axis to adjust the total angular momentum of the system. It is anticipated 
that the action of the nutation damper and momentum wheel under external torques 
will increase the total momentum of the system and necessitate periodic "dumping" 
by the MASC. 

The partitioning of the total momentum between wheel and body is achieved 
through the action of the pitch control electronics (PCE) in either a closed-loop 
mode (despun or 4 rpm) or open-loop mode (1 to 8 rpm). Details of command 
generation are given in another document. ^ 

2.4 CONTINGENCIES 

Two contingencies have been analyzed extensively: 

• Recovery from zero angular momentum failure. 

• Torquing to check out the Extreme Solar Ultraviolet Measurement 
(ESUM) field of view. 

^Computer Sciences Corporation Memorandum, Pitch Control Electronics 
Commands, R.D. Woolley, September 21, 1973. 


2-9 



2.4.1 Recovery From Zero Angular Momentum Failure Mode 


Prior to separation, the Delta launch vehicle is scheduled to spin up the AE 
spacecraft to about 10 rpm to impart the required angular momentum to stabi- 
lize the system. Failure to spinup (the zero angular momentum failure mode) 
results in serious consequences which endanger the mission. 

® With the spacecraft tumbling, the experiments cannot function 
® Attitude data s severely impaired 

• No orbit adjustment maneuvers can be performed 

It should be noted that the actual angular momentum will never be exactly zero, 
but will always have some finite value. 

The following procedure will be used to determine the angular momentum (mag- 
nitude and direction) if this failure to spinup occurs. 

• Spinup wheel to 10 radians per second 

• Determine from commanded wheel and measured body rates 
from Sun times and magnetometers 

• Determine nutation angle, 6 , from Sun angles and magnetometers 

• Determine total angular momentum from H = H /cos 0 . 

I u 

• Determine angular momentum direction with Sun and magnetometer 
data. 

Three recovery methods are possible: 

• Off-loading one tank with a delta-V thruster 

• Yaw thruster firing into a flat spin with passive nutation decay, 
using the momentum wheel to ensure positive spin 
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• Yaw thruster firing into a flat spin with active ground-controlled 
real-time nutation damping using the momentum wheel and mag- 
netic coil 

These methods are described in the following subsections. 

2.4. 1. 1 Off-Loading One Tank With the Delta- V Thruster 

If fuel is drawn continuously from the number 4 tank (along the -y-axis), the 
center of mass shift gives a moment arm along -y , providing a torque along 
the z-axis (Figure 2-1). 


T z » -Ay x F x 

Because the center of mass shift is linear in time, the angular momentum in- 
creases quadratic ally. Table 2-2 shows the values at 100-second intervals, 
where the impulsive approximation has been used for AH : 

AH = -Ay F At 

An additional source of angular momentum accrues from the momentum carried 
off by the fuel. Because the body can be made to spin in reverse by command- 
ing the tachometer loop on, subtracting off negative angular momentum will 
increase the total momentum of the system. 

AIT 2 
AH = -mr to 

z 

The values were calculated at the maximum tachometer rate of 41 radians per 
second. 
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It is estimated that approximately 1034 seconds of fuel usage (15.7 pounds of fuel) 
is required to add 318 in. -lb-sec of angular momentum, which will allow the 
spacecraft to be despun at the lowest tachometer rate of 10 radians per second. 

If time is critical, the burn should be performed as quickly as possible. If 
time is not critical, the burn can be split into three firings, and the following 
operational procedure is suggested: 

Determine tumble contingency from magnetometers and DSAI 
Command tachometer loop at 41 radians per second 
Load and fire 345-second continuous delta-V burn 
Determine attitude from magnetometers and DSAI 
Run MSAP for coil commands 

Load and execute coil commands to negative orbit normal 

[WHS and BHS will acquire the Earth. The nadir programmer will work at a 
negative spin rate. ] 

Load and fire second 34 5 -second delta-V burn 
Confirm attitude from WHS-Sun data 

[(DAPS becomes operational (pulsed mode) after nutation has decayed to less 
than 5 degrees to avoid the risk of the BHS missing the Earth and confusing 
the nadir programmer. 3 

Load and fire third 345-second delta-V burn 
Confirm attitude and rates 

[OAPS becomes operational in despun mode. If necessary, the orbit can be 
corrected using the opposite tank, and the angular momentum will increase 
further. 3 

Run MSAP for spin coil commands 
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Table 2-2. Angular Momentum From Off-Loading One Tank With Delta-V Thruster 


t 

F 

(POUND) 

F 

(POUND) 

Ay 

(INCH) 

Ay 

(INCH) 

AH j- 

rxF 

(IN.-LB-SEC) 

L O 

z 

(RAD/SEC) 

AH 2 
mr cj 2 

(IN.-LB-SEC) 

H TOTAL 
(IN.-LB-SEC) 

0 

3.983 


10' 14 






100 

3.916 

3.950 

-0.0156 

-0.0078 

3.081 

-1.3469 

4.774 

7.855 

200 

3.842 

3.879 

-0.0311 

-0.0234 

9.058 

-1.3388 

4.665 

21,578 

300 

3.750 

3.796 

-0.0463 

-0.0387 

14.690 

-1,3246 

4.529 

40.797 

400 

3.650 

3.700 

-0.0615 

-0.0539 

19.943 

-1.3048 

4.354 

65.094 

500 

3.551 

3.600 

-0.0765 

-0.0690 

24.843 

-1.2796 

4.156 

94.094 

600 

3.459 

3,505 

-0.0913 

-0.0839 

29.407 

-1.2497 

3.949 

127.450 

700 

3.376 

3.418 

-0.1061 

-0.0982 

33.730 

-1.2152 

3.741 

164.921 

800 

3,303 

3.340 

-0.1207 

-0.1134 

37.870 

-1.1765 

3.535 

206.295 

900 

3.238 

3.270 

-0.1352 

-0.1280 

41.846 

-1.1338 

3.333 

251.474 

1000 

3.181 

3.210 

-0.1496 

-0.1424 

45.703 

-1.0871 

3.132 

300.309 

1034 

3,165 

3.173 

-0.1544 

-0.1520 

16.657 

-1.0166 

0.976 

317.942 


[Leave attitude at a negative rate; run the spin coil to decrease this rate. ] 

Command tachometer loop at 10 radians per second 
Spinup to 1080 in. -lb-sec 

[Will require about 3 days at 24 in. -Ib-sec per orbit. 3 

The three firings can be accomplished at successive perigees to minimize the 
risk of changing the orbit. 

2. 4. 1,2 Yaw Thruster — Passive Damping 

This technique uses very little fuel (less than 0.5 pound) because of the 15-inch 
moment arm of the yaw thruster. The yaw thruster is fired for about 18 sec- 
onds to achieve slightly less than the nominally required angular momentum of 
1200 in. -lb-sec. Then the spacecraft is placed in 10 radians per second 
tachometer mode to maximize the rate of nutation damping toward the positive 
Z-axis. After the nutation half cone angle has decreased to less than 5 degrees 
(from the initial value of 90 degrees) and after the angular momentum vector 
has been maneuvered to the negative orbit normal using magnetics, orbit adjust 
maneuvers can be performed in either a spinning or a despun mode. 

If the perigee altitude is high enough so that aerodynamic pumping of nutation is 
small, the technique results in an operational spacecraft capable of orbit ad- 
justment maneuvers within 2 to 3 days. If the perigee altitude is below 150 
kilometers, then aerodynamic nutation pumping might significantly increase the 
time required for this technique, or perhaps even prevent the eventual success 
of this technique. The following procedure is required to implement this 
method: 

[pitch control electronics is off, so the wheel is not spinning significantly rela- 
tive to the body. 3 

Determine tumble contingency from MAGS and DSAI 
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[Verify that perigee is above 150 kilometers. ] 

Fire yaw thruster for 18 seconds to achieve approximately 1100 in. -lb-sec 
angular momentum 

Command spacecraft into tachometer loop at 10 radians per second 

Wait (1/2 to 1 day) for nutation to decay to about a 30-degree half cone angle 
[Determine attitude while waiting. ] 

Load and execute coil commands maneuvering to the negative orbit 
normal 

Wait for nutation to decay below a 5-degree half cone angle. 

Load and execute MASC commands to trim the angular momentum to 
1200 in. -lb-sec. 

2.4. 1. 3 Yaw Thruster — Active Damping 

This method could result in an operational spacecraft within 1 to 2 days, and 
uses less than 0.5 pound of fuel. However, it requires the development of 
nutation damping computer programs to operate in the Atmosphere Explorer 
Operations Control Center (AEOCC), phasing real-time commands about every 
5 seconds with quantities calculated from real-time telemetry. The efficiency 
of this technique in speeding up nutation damping over the passive damping case 
depends on the command timing accuracy (fractions of a second) achieved by 
that presently undeveloped software. 

After the nutation half cone angle has decreased to less than 5 degrees and after 
the angular momentum vector has been maneuvered to the negative orbit normal 
using magnetics, orbit adjust maneuvers can be performed in either a spinning 
or a despun mode. The time required before this operational capability is 
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possible is 1 to 2 days. The following procedure is required to implement this 
method: 

[Pitch control electronics is off, so the wheel is not spinning significantly rela- 
tive to the body. ] 

Determine tumble contingency from magnetometers and DSAI 

Fire yaw thruster for 18 seconds to achieve 1100 in. -lb-sec angular 
momentum 

Command spacecraft into tachometer loop at 10 radians per second 

Initiate active nutation damping maneuvers throughout every real-time 
station contact until the half cone angle is below 5 degrees 

As soon as the half nutation cone angle is below 30 degrees, load and 
execute (between station contacts) coil commands maneuvering to the 
negative orbit normal 

As soon as the half nutation cone angle is below 5 degrees, load and 
execute MASC commands to trim the angular momentum to 1200 in. -lb-sec 

2.4.2 Torquing to Check Out the ESUM Field of View 

The ESUM ejqjerimenter has requested that the spin axis be torqued through 
Sun angles of 35 to 90 degrees to check out the instrument's field of view. 

The maneuver is requested within 6 weeks of launch to be performed in about 
5-degree Sun-angle increments with the attitude determined and experimental 
data taken at each step. An analysis was made of magnetic maneuvers at two 
possible times: 

• first inversion from negative to positive orbit normal 

• six w'eeks from launch on an uncorrected orbit 
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2.4.2. 1 Inversion Time 


Although the time of first inversion on the current mission plan will occur 8 or 
9 weeks after launch, these plans may be revised or the experimenter's time 
requirement may be relaxed. Because a magnetic inversion maneuver is al- 
ready required at this time, torquing toward the Sun would not impose any 
additional operations nor any additional time outside mission constraints. 

Preliminary analysis established that maneuvers at declinations south of the Sun 
provided more WHS data than those north of the Sun. When torquing was for- 
bidden for 20 minutes on either side of perigee, the maneuver required 4 days, 
which seemed too long to have the spacecraft outside mission constraints. 
Therefore, continuous torquing was used to speed up the maneuver. 

To get a somewhat realistic orbit, elements were taken from an EPHEM tape 
for December 28, 1973 (from a November 1, 1973 launch) corrected for the 
difference in the node for the current December 7, 1973 launch date, giving an 
effective epoch at February 2, 1974. These elements were propagated by 
analytic routines (BRRS or FASTOG) to the maneuver time about February 13, 
1974. The perigee altitude is about 145 kilometers. 

For MSAP, nominal values were used for: 

• drag coefficient, C D = 2.0 

• angular momentum - 1200 in. -lb-sec 

• CP-CM offset = 0.1 inch along Z-axis 

• two coils, 98,000 pole-cm each 

• night coil reduction = 5 percent 

• initial attitude near the negative orbit normal 

Selected points (0 = 130 degrees, 90 degrees, 50 degrees) were input 

max 

to ODAP to check data availability and shadow times, which are shown on 
Figure 2-2. With the torquers on continuously, the Z magnetometer will be 
saturated, and it will be especially important to have the spacecraft in sunlight 
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RIGHT ASCENSION (DEGREES} SUN ANGLE (DEGREES) 



to get a second reference vector for accurate attitudes. The station coverage 
around the middle orbit (February 13, 1974 0437-0646) illustrates a potential 
problem in that the station coverage coincides with the times of missed data. 
Because attitude data might be severely degraded by blind orbits, it would be 
necessary to either turn on the tape recorders near apogee or to alter the atti- 
tude trajectory or maneuver times. 

2 . 4 . 2 . 2 Six Weeks After Launch 

Approximately 6 weeks after launch on an uncorrected orbit, the Sun angle will 
be at 35 degrees in a favorable geometry. A magnetic maneuver can be made 
away from the Sun to 90 degrees and back in about 36 hours. Figure 2-3 illus- 
trates the variation of Sun angle and right ascension with time. At this partic- 
ular geometry there is no shadow and no time when both WHS sensors miss the 
Earth (though WHS1 will miss at apogee and WHS2 will miss at perigee at the 
Sun angle of 90 degrees). The Sun angle is changing at about 7 degrees per 
orbit, but the torques can be arranged to have a different set of Sun angles at 
perigee on the two legs of the maneuver. 
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RIGHT ASCENSION {DEGREES) SUN ANGLE {DEGREES) 


90 - 
80 - 
70 - 
60 - 
50 - 
40 - 
30 - 


SUN ANGLE VERSUS TIME 
RIGHT ASCENSION VERSUS TIME 
NO SHADOW; NO MISSED DATA 
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Figure 2-3 


ESUM Torquing 6 Weeks After Launch 


SECTION 3 - ATTITUDE SUPPORT SYSTEM OVERVIEW 

3.1 MSAD EXECUTIVE AND GRAPHICS 

The Multi-Satellite Attitude Determination (MSAD) System has been developed 
for general attitude determination. This executive system employs an inter- 
active control structure to serve as the operational environment or base on 
which the AE-dependent attitude determination system was built. It is com- 
posed of a set of drivers and utility routines to control the sequence of appli- 
cation subsystems, graphic displays, core storage allocation, selection of 
routines, and error-handling capabilities. MSAD is used in conjunction with 
the MSAP System for complete attitude support. The total system is referred 
to as the Multi-Satellite Attitude Program System (MAPS). The MSAD/AE 
System is composed of a Simulator Subsystem, a Telemetry Processor Sub- 
system, and an Attitude Determination Subsystem. 

3.1.1 Simulator Subsystem Overview 

The Atmosphere Explorer Simulator Subsystem (AESIM) is designed to simulate 
the telemetry data required for attitude determination. During the prelaunch 
period, the primary function of AESIM will be to test the Telemetry Processor 
and the Attitude Determination Subsystems; during the postlaunch period, 

AESIM will be used as an analytical tool in the study of attitude dynamics. In 
addition, AESIM will create on demand a disk file of simulated data to be used 
in the prelaunch testing of the Operations Control Center (OCC) and the input 
processor (IP) data links. 

AESIM consists of two separate attitude simulators, TELSIM and ADSIM, each 
one tailored specifically for one of the functions mentioned above. TELSIM is 
more suitable in testing the telemetry processor and data links. It makes use 
of a simplified, non-nutating spacecraft model with either a constant attitude 
or a time varying attitude read from an MSAP tape. In addition, small-range 
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attitude maneuvers can easily be simulated with TELSIM by specifying the 
initial and final attitudes and the time span of the maneuvers. 

The Attitude Dynamics Simulator (AD SIM), on the other hand, simulates attitude 
data from a nutating satellite by numerical integration of Euler’s equations. Its 
inputs are the moments of inertia along the principal axes, the angular momen- 
tum vector, the body axes relative to the angular momentum vector, spacecraft 
position, and relevant data describing the spacecraft's sensor complement. 
ADSIM is also able to simulate coning, which is caused by an offset between 
the geometrical and principal axes of the spacecraft. 

3.1.2 Telemetry Processor Subsystem Overview 

The basic function of the telemetry processor is to convert raw data in telemetry 
counts from the AE-C spacecraft to geometric entities used in attitude determi- 
nation. During the conversion process, telemetry data resides in the telemetry 
processor in three forms: working arrays containing raw data in telemetry 
counts; intermediate arrays of sensor readings in engineering units; and output 
arrays of cone angles, dihedral angles, and azimuth determination information. 

The telemetry processor is called by TPDRIV, a driver built into the'MSAD 
executive system. Input to the telemetry processor consists of either real or 
simulated telemetry data contained on disk files. Additional input may include 
NAMELISTS, sensor data sets consisting of previously processed data in engi- 
neering units, and ephemeris information. Output consists of cone angles, 
dihedral angles, and azimuth determination information. Additional output 
optionally includes a data set containing accelerometer and OAPS data and a 
data set containing sensor information in engineering units. 

The telemetry processor processes data in both the real-time and definitive 
modes. The principal difference between the two types of processing is in the 
handling of output arrays. In the real-time mode, output arrays are passed to 
the Attitude Determination Subsystem when they are full. In the definitive mode, 
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full output arrays are written on a data set. Thus, in the real-time mode, 
several calls to the telemetry processor may be required to process a complete 
pass of data, whereas in the definitive mode only one call is required. 

The basic flow through the telemetry processor is to read a user- specified 
number of raw data records into a working array, unpack the data from the 
working arrays and store it in arrays of sensor readings in engineering units, 
then convert the data in the sensor arrays to cone angles and store these in 
output arrays. This cycle continues until an output array is filled; processing 
then proceeds as explained above, depending on whether the system data is in 
the real-time or the definitive mode. 

The telemetry processor subsystem is controlled by an executive which calls 
individual modules to perform specific functions. The executive also controls 
looping within the processor. Backward moves to the telemetry processor for 
reprocessing, either from another subsystem or from within the telemetry proc- 
essor, always return to the executive. Therefore, reprocessing starts at the 
beginning of an interval rather than in the middle of a partially processed inter- 
val. Reprocessing may start with raw telemetry or with sensor data in engi- 
neering units if a sensor data set has been written. 

Optional processing available in the telemetry processor includes quality byte 
checking, time checking, range checking, magnetometer bias determination, 
addition of constant biases to sensor readings, sifting of high- volume data, and 
display of data at various steps in the processing. 

The telemetry processor can run in either the interactive or the noninteractive 
mode (MSAD directory array mode). In the interactive mode, operator inter- 
vention is required only at the initial control parameter display to minimize 
attention handling in real-time processing. However, numerous optional dis- 
plays are available for analysis. These include plots and tables of all inter- 
mediate sensor arrays and all output arrays. The intermediate sensor array 
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displays include displays of all data points for a short period of time or selected 
data points for a longer period of time. There are also additional control pa- 
rameter displays which allow operator control over the optional processing 
features. 

In summary, the telemetry processor performs the required functions of read- 
ing raw telemetry, converting telemetry counts to engineering units, and con- 
verting these readings to cone and dihedral angles. Optionally, it performs 
various quality checks, displays data, and allows data to be changed or re- 
processed. 

3.1.3 Attitude Determination Subsystem 

The Attitude Determination Subsystem (ADS) has two major functions. In the 
real-time operational mode, it computes the attitude of the spacecraft's angular 
momentum vector and writes this information into the active file together with 
the averaged spin rates of the wheel and the body. This data set is used for 
control purposes. 

In the definitive mode, it computes certain attitude parameters on a uniform 
time-step basis and writes these parameters in a file for the experimenters. 
Using the data in this file, a three-axis attitude can be computed for any time 
T during the time span of the solution. 

The two operational modes are interchangeable on the same batch of data, giving 
either type of output. The mode cannot be changed, however, once the operation 
has started. 

The ADS consists of three modules: the filtering module (FILTER), the angular 
momentum module (ANGMOM), and the azimuth module (AZIM). 

3, 1. 3. 1 Filter Module 

The main functions of the filtering module are to reduce data volume, smooth 
data to remove random noise, and determine and remove the effects of certain 
types of dynamic motion from the data (precession, coning). 
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The filtering module processes the cone angle arrays passed from the telemetry 
processor. There are three types of these cone angle arrays: wheel horizon 
scanner arrays, Sun sensor arrays, and magnetometer arrays. 

In attitude computation, any combination of these sensor data can be used, and 
any data type can be skipped or downweighted at the user's option. Another op- 
tion is the sampling time interval for the filter output. 

These data arrays are processed by a recursive least-squares filter, which 
passes the parameters on to the next modules at equal time intervals. The 
processing of the data takes place point by point. The filter has a fading mem- 
ory, allowing any slow changes in the attitude parameters to be followed. The 
model function to describe the sensor cone angles has three parts: one term 
to describe the cone angle between the angular momentum of the spacecraft and 
the reference vector, a sinusoidal term to describe the motion of the Z-axis 
around the angular momentum vector {precession), and a second sinusoidal term 
for motions that show a periodic feature around the Z-axis with the body period 
(coning). This model function has seven parameters to fit the data. The ini- 
tialization of the filter is done by an application of the Fast Fourier Transform. 

The filter module also has simpler operational modes. It can be used for N 
point data averaging (running average mode) or data sifting if required. The 
recursive filtering mode is used mainly for precession determination in the 
definitive mode. 

The filter module also computes a dihedral angle array at the sampling times 
to permit resolution of attitude ambiguities in the angular momentum module. 

For each type of cone angle array there is a corresponding dihedral angle array. 
These angles are computed with respect to a fixed axis on the body so that the 
dihedral angle between two sensor events will be simply the difference of these 
dihedral angles. There is no adjustment for precession on these angles because 
the only purpose is to resolve ambiguity. The smoothed cone angle arrays and 
the raw dihedral angles are the data base of the angular momentum module. 
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3. 1. 3. 2 Angular Momentum Determination Module 

The capabilities of the angular momentum module are to determine the attitude 
of the angular momentum vector of the spacecraft; give one-point attitude in the 
real-time mode, which represents an average attitude for the pass or the up-to- 
date attitude after a perigee pass; provide an attitude history for the whole pass 
in the definitive mode; compute an initial attitude value — based on cone and di- 
hedral angles — to initialize the least-squares scheme (quick- look mode); batch 
or recursively process sensor data; and to compute third-axis attitude in the 
definitive mode. 

In the definitive mode, the preprocessed (filtered) data is written out in an 
auxiliary file if current definitive ephemeris data is not available. In the real- 
time mode, the filtered data array is processed immediately using the most 
recent ephemeris data. In the first step, the cone and dihedral angle data are 
used for an initial attitude determination; this is called the quick-look attitude 
determination mode. This can be done point by point or in batch mode. The 
output attitude is used to initialize either a batch (DCCONS) or recursive 
(GRECRS) least-squares differential correction subroutine to compute a refined 
attitude based entirely on the cone angle data. After the one-point attitude has 
been determined and written in the archive file together with average wheel 
and body rates, the operation of the real-time mode ends. 

In the definitive mode, as soon as ephemeris information is available, the fil- 
tered cone angle data goes through the angular momentum module to obtain 
angular momentum history. This procedure is identical with the one used in 
the real-time mode except that the attitude is computed at equal time steps and 
stored for further processing. 

The next step is the determination of the spacecraft Z-axis, The computation 
is based on the wheel horizon sensor data, with the Sun sensor data as an 
alternative. The subsystem is flexible enough, however, to use any available 
data type for Z-axis computation. 
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Even if both the coning and the precession components are present and separated 
from the data in the filter module, only the one having the larger contributing 
bias is considered for Z-axis computation. The Z-axis information is computed 
at equal time intervals, with about one point each body spin period. 

3. 1. 3. 3 Body X-Axis Determination Module 

In the definitive mode for three-axis attitude determination, the azimuth angle 
is also needed. This computation takes place in the azimuth module. The 
module's capabilities are to compute azimuth angles and rates at equal time 
intervals, and provide attitude parameters for experimenters, with up to 
quadratic interpolation terms. 

The primary source of data is the body horizon sensor readings in the spinning 
mode and the split-to-index pulses in the despun mode. In the despun mode, 
the pulses are converted directly to phase angles in the telemetry processor, 
and these angles are smoothed in the filter module. 

The Sun event times are used for calibration and the magnetometer data for 
backup. The attitude parameters are written out at equal time intervals, and 
a quadratic interpolation can be applied to compute three-axis attitude on a 
continuous basis. 

3. 2 MSAP SYSTEM OVERVIEW 

The AE version of MSAP is designed to generate the attitude control commands 
necessary to maintain spacecraft attitude within specified constraints. It is 
executed in the interactive graphics mode, with program flow under optional 
operator control. Computed spacecraft commands are both printed and written 
on a nonvolatile disk-file archive dedicated to AE for ease in subsequent han- 
dling. 

The MSAP/AE system has been designed to provide a flexible attitude prediction 
and optimal control system that may be expanded to support additional satellites. 
The system's essential features are described in the following subsections. 
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There are two archives used by MSAP/AE: Command History and Attitude 
History. 

The determined attitudes (along with appropriate times) are in the Attitude 
History archive and are entered by the Attitude Determination Subsystem, pro- 
viding the link between the attitude determination and control functions. 

The Command History archive contains sequential entries of two types: 

1. spacecraft commands affecting attitude control 

2. determined values of changing parameters (x, y, z CP-CM offsets) 
essential to computing required commands 

These two archives are read by setting two FORTRAN NAMELIST parameters 
enabling the operator to streamline initialization for parameter estimation. 

Spacecraft attitude control commands (and their times) are determined by 
MSAP/AE and written into the archive for use in command displays and in 
subsequent simulations with overlapping time periods. Canceled commands 
are not deleted from the archive; their canceled status is indicated by entry 
of a flag which allows reinitialization of their initial types and times for con- 
trol center use. 

Estimated values of the center-of-mass location, moment of inertia, and 
expected changes of these are included in the archive. 

The archive is allowed to accumulate entries until near full. A utility program 
periodically erases old entries. 

3.2.1 Input Subsystem 

All input to the system is accomplished using the FORTRAN NAMELIST feature. 
In addition to the FORTRAN NAMELIST, the history data and the observed 
attitude data will be inputted as formatted tables which can be filled from the 
Command History and/or Attitude History archives. 
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3.2.2 Attitude Prediction Subsystem 


The attitude prediction subsystem simulates changes in spin-axis attitudes and 
spin rates. The simulation algorithms are based on the assumption that nuta- 
tion is of negligible magnitude, and that collinearity is therefore maintained 
between the spin-axis direction and the spacecraft angular momentum vector. 
Attitude-dependent disturbance torques are computed in the simulation and 
numerically integrated to predict changes in the spacecraft angular momentum 
magnitude and direction. 

All disturbance torque calculations are made at user option. Included are grav- 
ity gradient, aerodynamic, magnetic, direct solar radiation pressure, and 
thruster torques. 

The predictor also monitors day /night conditions to determine whether or not 
there is solar radiation pressure to be computed and to determine the magnetic 
moment strength during low voltage conditions. 

It also simulates momentum transfers between the main body and wheel, as it 
is controlled by the onboard pitch feedback-control system. 

3.2.3 Attitude Control Subsystem 

This subsystem automatically computes the optimum commands subject to 
operator-specified constraints to precisely implement control functions. The 
types of commands automatically computed by this subsystem are as follows: 

• time-optimal or energy-optimal spin-axis coil commands to bring 
the spin axis to a specified attitude 

• momentum coil commands to bring the angular momentum magni- 
tude to a specified value 

• yaw thruster commands to bring the spin axis to a positive or neg- 
ative orbit normal 
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• auto-roll commands to maintain the attitude near the orbit normal 

• tachometer commands to achieve the desired body spin rate 

3.2.4 Output Subsystem 

The output subsystem writes the set of computed commands onto the AE Com- 
mand History archive, produces printed output, and optionally writes an attitude 
history tape. 

3.2.5 Parameter Estimation Subsystem 

The parameter estimation subsystem is designed to provide a means of estimat- 
ing nonmeasurable parameters necessary for torque computation, such as the 
value of the CP-CM offset, or the value of a residual magnetic dipole moment, 
etc. Such values are necessary to generate valid attitude predictions using the 
attitude prediction subsystem and are consequently necessary to plan a com- 
mand sequence correctly. 

This subsystem, KOPTIM, uses an operator-selected sequence of attitude de- 
terminations and spacecraft commands from the MSAP archives, then itera- 
tively determines the values of a set of operator-specified MSAP parameters 
which minimize the weighted squared error between a set of observations and 
predictions at corresponding times. KOPTIM performs the following functions: 

1. KOPTIM inputs information from MSAP that specifies which param- 
eters are to be optimized. Then it converts these parameters into 
the state vector from subroutine KOPTM1, which is a pure function 
optimization program. 

2. KPTIM calls subroutine KOPTM1, which generates a new set of 
state vectors according to the least- squares method. 

3. KOPTIM converts the new set of state vectors to the parameter 
values of MSAP. 
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4. If the program, satisfies the convergence criterion or exceeds an 
iteration limit, KOPTIM sends a signal to MSAP to terminate the 
computation. 

The optimization technique used by KOPTM1 is a first-order adaptive method 
to find a local minimum of the nonlinear objective function. The original 
method was developed by Gauss and Newton and is obtained from a Taylor 
series expansion of the model through the linear terms. Using the idea of 
Levenberg to control convergence , Marquardt and Kotsakis have introduced 
modifications to the basic method to achieve a nearly optimum convergence 
performance. 

The advantages of KOPTM1 are as follows: 

1. The objective function decreases monotonically from iteration to 
iteration. This implies that the set of parameters computed last 
is a better set over the previous iterations. This is an advantage 
if the computation has to be terminated prior to the ultimate con- 
vergence. 

2. An optimum interpolation between the gradient and Gauss-Newton 
method is chosen adaptively by the algorithm whereby a proper 
direction and magnitude for the correct vector is automatically 
computed to yield rapid convergence. (Mathematically this is 
controlled by the value of A. in Equation (3-3). ) 

A brief summary of the KOPTM1 algorithm is as follows: 

Let the prediction vector to be fitted to a set of data points (t. Z.) , 

j = 1, 2, . . . , n be 

Y = f(t,X 1 ,X 2 , ... ,X np )-f(t,X) (3-1) 
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where X is a parameter vector to be estimated. Also, let the weighted least- 
squared function be 


Q (X) = [Z - Y] T W[Z - Y] 


(3-2) 


The KOPTM1 algorithm to minimize this objective function Q (X ) is the re- 
cursion 




6x. = fj T WJ. + x.l I " 1 j T w[Z - f (t , X.)] 
-iLi i iJ i — - — l 


(3-3) 




i = 1, 2, . . . , np 


where J. is the nxnp Jacobian matrix of the first partial derivative of 
1 f(t , X) evaluated at X. and X. is a scalar, A. ^ 0 . 


It is shown that any value of A. such that 


Q (x. + 5x (\)) < Q (X.) 


(3-4) 


will ensure convergence and the choice of X. will affect the rate of convergence. 
Various strategies have been proposed to determine an ''optimum’ 1 value of A 
dynamically at each iteration. 

2 

KPTM1 uses the scheme modified by Marquardt and Kotsakis to control the 
convergence of Equation (3-3). 


2 

Computer Sciences Corporation, 9101-12300-02TM, Data Management in 
Atmosphere Explorer Support Software , G. Kotsakis, January 1973. 
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The convergence criteria used in KOPTM1 are as follows: 


xf> -xf- 1 * 


x. <K) -x! K ' 1) l 
1 1 


(K-l) 


<P 


<p' 


for i = 1, 2 


(3-5) 


for i = 3, 4, . . . , np 


where K = iteration index 

X , X = initial values of right ascension and declination (radians) 

X 2 

p , p 1 = the order of accuracies desired for convergence 

The maximum number of KOPTM1 iterations depends on the accuracy desired 
for the convergence. However, five iterations should be sufficient for most 
cases. One KOPTM1 iteration requires (np + 2) function evaluations (MSAP 
predictions) , where np is the number of parameters being optimized. 
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SECTION 4 - OPERATIONS PLAN 


This section presents nominal procedures and schedules for several operations: 

• launch sequence 

• perigee lowering 

• apogee restoration 

• normal operations 

• low perigee operations 

• magnetic inversion 

• yaw thruster firing 

• system momentum determination 

4 . 1 LEAD TIME GROUND RULES 

The following ground rules have been established for the lead time required for 
commands being sent through AEOCC: 

• In general, a 6-hour lead time is required 

• A 10- to 12-hour lead time is required before each OAPS bum 

• Only a 1.5-hour lead time is required for apogee restorations when 
performed sequentially (on every second or third perigee) 

On the schedules presented on the following pages, station coverage is indicated 
by orbit number, which is incremented each time the spacecraft goes through 
the ascending node. The letter P indicates perigee passes and A indicates 
apogee passes. The orbit number is indicated by the numbers following the P 
and the A. These numbers are representative. 
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NOMINAL LAUNCH SEQUENCE 


STATION 

SPACECRAFT EVENTS COVERAGE 


ATTITUDE OPERATIONS 


ATTITUDE DETERMINATION (AD) USING DSAI-BHS-MAGS; 

CHECK SPIN RATE USING MAGS 

AD; DETERMINE NUTATION AMPLITUDE 

REPROCESS WITH UPDATED ELEMENTS; ATTITUDE PREDICTION (AP) 
FOR QOMAC AND MASC COMMANDS TO BE LOADED ON A2 


LAUNCH, INJECTION, YAWTO NEGATIVE 
ORBIT NORMAL SPINUP TO 10 RPM 

SEPARATION 

(SANTIAGO PASS) I PQ — \ 


PRELIMINARY ELEMENTS 

(TANANARIVE PASS) 

AO 1 

UPDATED ELEMENTS 

) 

I 

1 A1 J 

unoi i i Mr t i 

SPINUP MWA IN TACH LOOP; LOAD QOMAC AND } 
liiaer rnMMflMnc 1 

PRECESS TO ORBIT NORMAL 

A2 | 


1 A3 1 

1 

1 

ci hsf pitch i nnp | 

1 M I 


LOAD QOMAC AND AUTOROLL COMMANDS | 

COMMAND MEMORY LOAD AND DUMP ^ 

1 A5 1 

1 A6 1 


AD; CHECK NUTATION DAMPIN G; DETERMINE MOMENTUM 
MAGNITUDE 

AD USING WHS-DSAI-MAGS; DETERMINE MOMENTS OF INERTIA 


AD; EVALUATE DRIFT RATES AND OA BIASES 

|“ AP FOR QOMAC AND AUTOROLL FOR LOADING ON A5 
AD; CHECK 4-RPM MODE 

AD; CHECK AUTOROLL PROGRAMMER 

AD; READOUT FOR ATTITUDE CONTROL 


COMMAND DESPIN 

LOAD QOMAC COMMANDS _ 

PRECESS TO OR BIT NORMAL 
OAPS CALIBRATION BURN 
PRECESS TO ORBIT NORMAL 

ORBIT TAPE 
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NOMINAL PERIGEE LOWERING 


STATION 

SPACECRAFT EVENTS COVERAGE ATTITUDE OPERATIONS 


LOAD QOMAC 


PREROLL TORQUES 

DESPIN 

FIRE OAPS 

AUTOROLL MEASURE 

SPINUP 

POSTROLL CONDITIONAL TORQUES 


A16 


A17 


| A18 | 
f A19 | 

' 

I A2 ° 1 

| A21 J 


AD 

AP FOR QOMAC FOR PREROLL AND AUTO- 
ROLL FOR POSTROLL FOR LOADING ON A19 


AD 


AD USING PREDICTED ELEMENTS; 
MESA/OAPS CHECK AUTOROLL 


AD USING UPDATED ORBITAL ELEMENTS; 
MOMENTUM DETERMINATION; AFTER 
SPINUP 
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NOMINAL APOGEE RESTORATION OPERATIONS (SEQUENTIAL) 


STATION 

SPACECRAFT EVENTS COVERAGE ATTITUDE OPERATIONS 


LOAD QOMAC 

PREROLL TORQUES 

DESPIN 

FIRE OAPS 

SPINUP 



A30 


\ 


AD 

AP FOR QOMAC PREROLL COM MANDS FOR 
LOADING ON A29 


AD 


AD 


AD USING ORBITAL ELEMENTS, AP QOMAC 
TO PREROLL FOR A32 


LOAD QOMAC 

DESPIN 

PREROLL TORQUES 
FIRE OAPS 

SPINUP 



A33 


AD USING EPHEMERIS DATA 


AD USING ORBITAL ELEMENTS; MOMENTUM 
DETERMINATION; AFTER SPINUP 
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NORMAL OPERATIONS 


STATION . 

SPACECRAFT EVENTS COVERAGE ATTITUDE OPERATIONS 


LOAD COI L COMMANDS 

COIL TORQUES TO 

OVERSHOOT ORBIT NORMAL 



1 A48 


A49 AD 


COIL TORQUES IF REQUIRED FOR 

PERIGEE BELOW 140 KM ( 

I A5 ° I 


I A51 I 

RECEIVE NEW EPHEMERIS 
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NOMINAL LOW PERIGEE OPERATIONS 


STATION 

SPACECRAFT EVENTS COVERAGE ATTITUDE OPERATIONS 


SPECIAL QOMAC LOAD 
PREROLL 


PREROLL 

LOAD QOMAC 


PREROLL 


AUTOROLL MEASURE 
AUTOROLL TORQUE 

PREROLL 


PREROLL — 

PREROLL 

LOAD QOMAC 

PREROLL 


PREROLL 

AUTOROLL MEASURE 

AUTOROLL TORQUE 
PREROLL 


PREROLL 


[ A41 | AD 

I AP FOR QOMAC COMMANDS FOR LOADING 

( ON A42 AND A44 



AD 



AD 


A44 


| P45 
| A45 


A46 


AD 

AP FOR QOMAC COMMANDS FOR LOADING 
ON A48 

AD 



AD 

AD 

AD 

AP FOR QOMAC COMMANDS FOR LOADING 
ON A52 

AD 

AD 
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NOMINAL MAGNETIC INVERSION OPERATIONS 


STATION 

SPACECRAFT EVENTS COVERAGE ATTITUDE OPERATIONS 


NEW EPHEMERIS 


LOAD QOMAC 


r 


CONTINUOUS 
COIL COMMANDS 
FOR 

APPROXIMATELY 
48 HOURS 


NEWEPHEM 


LOAD QOMAC 


NEW EPHEM 


V 


A61 | AD; CHECK ANALYSIS FOR MANEUVER TRAJECTORY 

f— AP; COIL COMMANDS CONTINUOUS FOR 24 HOURS 


I I 

| A63 | 
| A64 | 
A65 | 
A66 


AD 


1 

1 

_( 

A67 1 

• 

• 

• 

f 

A ? 2 | 

* ) 

• [ 

• J 


r a75_ i 


AD 


• J — AP COIL COMMANDS FOR 24 HOURS 




A79 AD 


A82 | AD 

j — AP FOR FINAL TRIM, IF NECESSARY 


A83 


I A84 I 


A85 


AD 
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NOMINAL YAW THRUSTER FIRING 

STATION 

SPACECRAFT EVENTS COVERAGE ATTITUDE OPERATIONS 



AD; DETERMINE CONTINGENCY CONDITION 
AP; THRUSTER COMMANDS 


AD 

— AP QOMAC COMMANDS TO ORBIT NORMAL 
AD 


AD 
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MOMENTS-OF-INERTIA DETERMINATION (AFTER PERIGEE LOWERING) 


STATION 

SPACECRAFT EVENTS COVERAGE ATTITUDE OPERATIONS 


LOAD DESPIN, OAPS, AND SPINUP COMMANDS — 
DESPIN — 

FIRE OAPS 
SPINUP — 


{ A 150 

— AD (USING PREDICTED ELEMENTS) 

— MOMENTUM DETERMINATION FROM WHEEL RATE 

1— MOMENTS-OF-INERTJA DETERMINATION FROM 

WHEEL AND BODY RATES 
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4.2 DETERMINING SYSTEM MOMENTUM 


4.2.1 Initial Determination 

Current launch plans call for the PCE to be off for the first couple of orbits. 
The wheel is expected to spin up to the body rate of about 10 rpm: 


H = Iw 


After nutation has damped down, most of the momentum will be about the 
z-axis: 


H TOT H 3 *33 ^3 


where I is well known and includes the flywheel, and oo can be determined 

oo o 

from Sun times (inertial), BHS leading edge times (Earth referenced), or 
magnetometer azimuth angles (B-field referenced). 

If there is a 1-percent error in I , however, there will be a 1-percent error 
in H tqt . When a subsequent despin is performed, this error will limit the 
accuracy of measurement of the wheel moment of inertia. 

4.2.2 Subsequent Determinations 

After a firing, the body moments of inertia will change and an accurate determi- 
nation of system momentum will require more than just wheel and body rates. 
The recommended procedure is as follows: 

• Despin the spacecraft 

• Determine the wheel spin rate 

• Calculate the total momentum 

• Command body to a high spin rate (4 to 6 rpm) 

• Calculate body (I ) moment of inertia from total momentum 
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SECTION 5 - DYNAMICS MODELING AND DATA CHARACTERISTICS 


5. 1 INVESTIGATION OF BODY DYNAMICS 

For a better understanding of the complicated dynamics expected on AE, the 
Continuous System Modeling Program (CSMP) was used to simulate the space- 
craft's motion. This package, which is available on the S/360-95, is especially 
useful because of the ease of development, modification, and output. The 
features listed below have been implemented. 

AE-C is modeled as a rigid body with 

• Momentum wheel, which can be offset from the Z-axis 

• PCE tachometer loop at a commanded rate, including electronic 
gains and delays 

• Fluid loop nutation damper, modeled as 10 concentric flow rings 

• External torques (only aerodynamic torques are computed at present) 

Aerodynamic drag force is computed inertially and rotated into body coordinates 
for torque computation: 

• Euler summetric parameters are integrated to provide the rotation 
matrix 

• Parameters are renormalized to improve accuracy 

• Drag force model has a constant inertial wind direction, an ampli- 
tude which is Gaussian in time, scaled according to the instantaneous 
area presented to the wind 

• Spacecraft is modeled as a cylinder with a variable center of pres- 
sure (CP) offset according to CP = 0.25 inch ±0.1 inch per 
degree yaw 

• Accelerations are computed at each integration step 


5-1 



To improve computational speed, there are branches around the computations 
of aerodynamic torque and accelerometer output when not required. 

Sample MSAD plots are shown of the body motion in Figures 5-1 through 5-4 
before the aerodynamic torque is applied. Figure 5-1 shows the cone angle 
between the angular momentum vector and the Z-axis versus time, Figures 5-2, 
5-3, and 5-4 show the Euler angles, 0 , <p , ij) , versus time. 

5.2 DATA CHARACTERISTICS 

A telemetry simulator (TELSIM) was built to test the Telemetry Processor and 
Attitude Determination Subsystem. It has the capability of simulating either 
telemetry frames or cone angles derived from the data. It can simulate various 
types of noise and bias, including the quantization of the telemetered values, 
which has been extremely useful in anticipating the actual data characteristics. 

Figures 5-5 through 5-9 are MSAD plots of typical expected data. Figure 5-5 
is a plot of the raw magnetometer cone angles before fitting (notice horizontal 
predicted line) of data. Figure 5-6 is a plot of the raw magnetometer cone 
angles after being smoothed. In addition, this plot show's the large change in 
cone angle at apogee caused by a bucket change in the Z -magnetometer read- 
ing. Note that Figures 5-5 and 5-6 are for two different time periods. Fig- 
ure 5-7 shows typical examples of a continuous batch of WHS data after being 
smoothed. Figure 5-8 is used to determine if the spacecraft is spinning in 
the proper direction. Figure 5-9 shows the single- versus the dual-scanner 
nadir angle solutions. It is possible to identify each scanner solution (plot) 
by operator intervention. 
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4 > (DEGREES) 





comrnzx m^cao 


DRTPL3 RAW VS. SMOOTHED MAG CONE ANGLES 



31.40.000 33.20.000 35.00.000 36.40.000 

73. 12. 14. 06 - 73. 12. 14.06 

XI TIMP X2 MRGT 

Y1 PRED T2 MAGO 

x= xf= r= tf= 


CPO I NT = 0RTPL3 WHAT NOW ? CALL OISPLAT D1SP 1 OF 1 


Figure 5-5. Raw Versus Smoothed Magnetic Cone Angles Before Smoothing 
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u^mrozD mzQo 


QATPL3 


RAW VS. SMOOTHED MRG CONE ANGLES 


1 DENT 


118.602; 
118.029- 
117.456 
116.883 
US. 310 
115.737 
115.165 
114.592 
114.019 
113.446 
1 12.873 
112.300 
111.727 
111.154 
1 10.581 


4 - 4 - 44 - 4 - 44 - 4 - 



4 - 4 - 4 - 4 - 4 - 


EXPAND 

REMOVE 

REPLACE 

REGEN 

RESTART 

FETCH 

STORE 

EXTFLAG 

INTFLAG 

EXTINT 

NO FLAG 


END 


1 10 . 008 J — * 

42.30.000 44.10.000 45.50.000 47.30.000 49.10.000 

43.20.000 45.00.000 46.40.000 48.20.000 50.00.000 


73. 12. 14.06 - 73. 12. 14.06 


XI TIMP X2 MAGT 

Ti PRED T 2 MAGD 

X= XF= T = TF = 


CPQ I NT =DRTPL3 WHAT NOW ? CALL DISPLAY OISP 1 OF 1 


Figure 5-6, Raw Versus Smoothed Magnetic Cone Angles After Smoothing 
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ojmrozD m^Qo 


94.721 

94,710 

94.698 

94,687 

94.676 

94.664 

94.653 

94.641 

94.630 

94.618 

94.607 

94.595 

94.584 

94. 572 
30.50.000 32.30. 00 

31.40.000 3: 

73. 12. 14.06 



XI TIMP 
Ti PRED 
X = 


X2 WHST 
T2 WHSD 



uorn^cn :z::d m^on 









cnnrn jocnmo— mrc nzi ji-odz 


TRIPLT SINGLE VS OUfiL SCANNER NADIR ANGLE SOLUTIONS 



73.12.1*1.06 - 73 • 12 . 1*1 . 06 

XI DNT X2 GNT X3 GNT 

T 1 GNAD1 12 GNAD2 T3 GNAD12 

X= XF= Y - YF = 


CPOI NT =NAD01 S WHAT NOW ? CALL DISPLAY DiSP 1 OF 2 


Figure 5-9. Single- Versus Dual-Scanner Nadir Angle Solutions 
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APPENDIX A - ORBITAL ELEMENTS FOR AE-C 


Epoch 

Semimajor axis 

Eccentricity 

Inclination 

Mean anomaly 

Argument of perigee 

Right ascension of ascending 
node 

Height of perigee 
Height of apogee 


06 hours 23 minutes GMT on 
December 14, 1973 

8614. 0 kilometers 

0.2414 

68.1 degrees ±0.1 degree {3a) 

359. 2 degrees 
166. 0 degrees 

251.4 degrees 

156 kilometers ±16 kilometers (3 a) 
4300 kilometers ±30 kilometers (3a) 
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APPENDIX B - AERODYNAMIC DRIFT 


The following printouts are representative of the predicted AE-C spacecraft 
drift behavior due to aerodynamic torques at a low perigee of 120 kilometers. 

The CM-CP offset is the radius vector from the center of mass to the center of 
pressure. This version of MSAP approximates the AE configuration as a right 
circular cylinder with the center of pressure at the centroid of the projected 
area, which is taken as the origin of the spacecraft reference system. The 
CM-CP offset in this case is due to having the CM location offset from the cen- 
ter of the reference system. 

Figure B-l shows a plot of the AE-C spacecraft declination versus right 
ascension due to aerodynamic torques with a CM-CP offset of 0. 1 inches at a 
low perigee of 120 kilometers. 
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DECLINATION (DEGREES) 






319 320 321 


RIGHT ASCE 


Figure B-l. 


Plot for One Day of AE- 
Ascension Due to Aeroch 



NSION (DEGREES) 


Spacecraft Declination Versus Right 
amic Torques 
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M S A P 


SATELLITE 


INDEPENDENT 


INPUTS 


KSTART 

• 

73 

12 

10 

5 S3 0 

PREDICTOR START TIME (YR » « MO. *DA. .HR..MIN.* SEC. 1 • 


KSTOP 

- 

73 

12 

1 1 

5 S3 0 

PREDICTOR STOP TINE < YR . .MO. .DA * .HR . .MI N. .SEC. > * 


MJMJNT 

* 




30 

NUMBER OF INTEGRATIONS PER OPBIT. 


KUNOUT. 

3 




30 

NUMBER OF OUTPUT POINTS PER CR0IT. 


MPERCC 

3 




132 

ORBITAL PERIOD IK MINUTES. 


MOTOR Q 

2 




90 

TORQUE COMPUTATION STEP SIZE <SEC*>. 


ALPHA 

3 




332.40 

INITIAL RIGHT ASCENSION OF SPIN AXIS (DEC.). 


DELTA 

-» 




-21.90 

INITIAL DECLINATION OF SPIN AXIS < DEG • ) • 


CMEGA 

= 




o 

o 

• 

<* 

INITIAL SPIN RATE OF SPIN AXIS <RPM.». 


RM AX 

S 




300.00 

ALTITUDE ABOVE MHICH AERODYNAMIC TOPCUE MAY BE NEGLECTED 

( KM* ) , 

GRACE 





0.0 

GROOMER DRAG TERM < RAC . / SEC . * *2 > * 


IPRINT 

m 




6 

STANDARD FPINTEO OUTPUT CSRN. 


INCARC 

= 




73 

CARD INPUT DSRN. 


ITAPE 

m 




10 

RESERVED FOR FUTURE OSRN. 


icpei 

= 




0 

OR B 1 TAPE DSRN. 


ioump 

s 




76 

INTERMEDIATE PRINTED OUTPUT DSRN. 


1EPHEM 

s 




0 

EPHEMERIS TAPE OSRN, 


MSG IN 





13 

ERROR MESSAGE DATA SET INPUT OSRN. 


msgolt 

* 




6 

ERROR MESSAGE OATA SET OUTPUT DSRN. 


CEPOCH 





7 312C7.00 

EPOCH DATE OF ORBITAL ELEMENTS (VYMMCD.Ol* 


CTIME 

- 




21160.00 

EPOCH TIME OF ORBITAL ELEMENTS I SECONDS > . 


GAX I S 

* 




EC25.12 

SEMI-MAJOR AXIS IKM.I. 


CEPSLN 





0.2421S3 

ECCENTRICITY. 


CIOTA 





€8.1 1 

1NCLINATICN CDEG . > • 


QUPSLN 

* 




359.41 

MEAN ANOMALY (DEG.). 


OONEGA 

3 




165.75 

ARGUMENT CF PERIGEE CCEG.l. 


C ALPH A 

* 




236.90 

RIGHT ASCENSION CF ASCENDING NODE (DEC.). 


OMOREF 

3 




731210.00 

REFERENCE DATE FCR LONGITUDE CF GREENWICH ( YYMMDD.O) • 



LCNGITUDE OF GREENWICH REFERENCED TO PREDICTOR START DATE (DEG.). 


OLAMCA 


78*53 
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I AT ARC 
I MS AflC 
I OFF < 1 i 
IOFF ( 2 ) 

I OFF <31 
I OFF ( A ) 

ICFF<5> 

< I FREE 

ISP*00 
IPANGL 
NSTEP1 
NST6P 2 
¥A TAPE 
IPERIN 

IPER 
HIKPEP 
*C F CFF 
YCMCFF 
2CMCFF 
XB IAS 
YB IAS 

2E FAS 

| 

WMH— - 
*M«-I** — ■ 


VHOXhS 



923203564*000 


BODY V MOMENT CF INERTIA (GM-CF**2) 
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ZMQINS 

1008275456.000 

BODY Z MOMENT OF INERTIA <GM-C***2! 

AMASSN 
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BODY MASS (GM) 

COIL1 
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COIL2 
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THRUST OF CELT A — V 1 ROCKET ILBS> 

OCLV2 

— 5-i rCW3- "■ 

— THPUST-OF PELTA-Vg ROCKET 
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GRAPH I C9 T TME" INC RE PENT C MINI 1 — 

ALT 
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IPSCDE 


PARAMETER OPTIMIZATION CODE 

CD 

2.000 

AERO DRAG COEFFICIENT 
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3.9071 

1196.1 

60.3 

167.7 

170,4 

e : 45 : 1 5 

e : 49 :4G 

323.3 

333.3 - - - 

-21.7 

-24,7 

3. s e?3 
3.-9-080 

1109.5 

1439,6 

60.7 

-6-0,7 - - • 

167.3 

467-^2 

179.0 
^ ^ 0^0 

e:s4: 5 

3 23.3 

-21 .7 

3.9e<n 

1189.6 

60.7 

166.9 

179.0 

e :5s: 3C 

323.3 

-21,7 

-3,9902 

1489.6 

-.60,7 

166.6 

170.0 

9: 2:55 

323.3 

-21.7 

3.9913 

1189.6 

60.6 

166.5 

179.0 

o: 7:arc 

323.3 

-21,7 - 

3.9922 

1189.6 

60.6 ... 

-166,5 

179,0 

9:11 :45 

323.3 

-21.7 

3.9929 

1109.6 

60.6 

166.6 

179.0 

— 9 : 1 6 :io 

3^-3, 3 - - 



a 1^ 

44^9,^ 

60 « 6 

16 6.7 

1 79 ■ 0 

9:20:3* 

323.3 

-21.7 

3.9440 

1 189. 6 

60.6 

166 .9 

179.0 


2 23.3 

-21,7 

.3,9 944 

1.13 9. ,6 

60,6 — 

_ 167.1 

-1X9*0 

9-5 29:25 

323.3 

-21.7 

3.9947 

1189.6 

60.6 

167.2 

179.0 



B-7 


Afr - C ATT t TUB e S U WM A WY- 


* frIGHT 
ASCENSION 

/ nprnFF^ i 
l U C vKUUo v 


DECLINATION 
(DEG R EES > 


SPIN RATE 
- 


MOMENTUM 


GAMMA 

- (peg , i 


LAMOA 

(OCG . l 


PHI MAX 
IO E G .I — 


10 73 
OEC 1C 73 
OEC 10 73 
OEC 10 73 
-DEC *0 73 
OEC 1C 73 
“DEC 10 73 
DEC 10 73 
-OEC 10-73 
DEC 10 73 
-Dee 10 73 
DEC 10 73 
OEC 10 72 
DEC 1C 73 
OEC 10 72 


9i 33 rfrfr 

9:38: is 
-■ 9:4*440 

9:47: S 
- - 9151 ! 30 
9:55:55 
10 : 0:20 
10: 4:45 
ID! 9210 
10:13:35 

10 : 16 : 0 
10 : 22:35 
10 : 26:50 
10 : 31 : 15 
I 0 : 35 : 4 c 


3 £3 *3 

223.3 

21 ^ 

-21.7 

: 3.9949- — - 

3.9950 

1 1 Ov .6' 
1 189.6 

frfr*6 

60.6 

167.4 

167.5 

— 179,6 

179.0 
179,6 
179.0 

223*3 

-21.7 

2i 9 95 1 
3.9551 

1199.6 

60.6 

60.6 

167.7 

m bJ ♦ J 

223.3 

-21.7 

3*9550 

1109.6 

60.6 

167.6 

1*79,0 — 
179.0 

2 23*3 — 

-21x7 

3*9949 - 

1-1 69.6- 

- 60.6 

W7.E- - 

179.0 

223.3 

-21.7 

3.5547 

1109.6 

60.6 

167.7 

179.0 

223.3 

—21*7 

-3.9 944 

1189.6 

60.6- - — 

l-67-.fr — 

179*0 

223.3 

— 21*7 

3.9541 

1 1 69.6 

60.6 

167.5 

179.0 

22 3.3- — 

- — 21V? 

— — -3.5536 

- -1189x6 

4 0.6 - - ■ 

167,4 

—179*6 

223.3 

-21 .7 

3.9931 

1189.6 

60.6 

167.3 

179.0 


2 23*3 
223.3 
223*3 


-24 *7 
-21.7 

-et •? 


3.9923- 

3.9515 

3*5504 


1169.6 

1169.6 

1189.6 


60.6 

60.6 

60.6 


-167^2- 
167.1 
— - 1 6 7 »-g - 


179.0 
176.9 
— • I70,fr- 


DEC 

10 

73 

10:40: 5 

223.3 

-21 . 7 

3.9693 

1 139.6 

60.6 

167.4 

178.9 

OEC 

10 

73 - 

10:44:30 

523.3 

at, 7 - — 

- 3.9882 

-4180.6 

60.6 

lfr7-.8 

— 1 70.9 

dec 

10 

73 

1 0:48:55 

323.2 

-21.7 

3.9674 

1109.6 

60.6 

168.3 

178.9 

DEC 

10 

73 

10 : 53 : 2 c 

323.5 

— 21 .6 

3.9870 - - 

1189.9 

69.7 

4fr7. 4 

- 1-78.7 

DEC 

10 

73 

10 : 57:45 

2 23.7 

-21.5 

3.9673 

1184.9 

61.0 

167.4 

178.5 

OEC 

10 

72 

11 : 2:10 

323.7 

-21.5- - 

-a. 48 91 

1185*0 

61.0 

— *• 467,3 

178.5 

DEC 

1C 

73 

11 : e:35 

223.7 

-21.5 

3.9692 

1185.0 

61.0 

167.1 

178.5 

oec 

10 

73 

11 : 11 : 0 

223.7- 

-21 .5 - 

3*9 90 3 

M-8fr*0- 

— 61.0 

.--■-407*0 — 

- 4 70*5- 

OEC 

10 

73 

11 : 15:25 

323.7 

-21.5 

3.9914 

1185.0 

61.0 

166.9 

178.5 

DEC 

10 

73 

11 :i9:so 

323.7 - 

-24.5 •- 

3.9923- — - 

-1105.0- 

64,0 

— 466.8 - - 

176.4 

OEC 

10 

73 

11 : 24:15 

323.7 

-21 . 5 

3.9530 

1185.0 

61.0 

167.0 

178.4 

OEC 

10 

73 

11 :2a:4c 

223.7 

-21.5 - - 

3.9936 

1 1 06.0 — - 

“61.0 

— 167.1 

-470.4 

DEC 

10 

73 

11 : 33 : 5 

2 23.7 

-21.5 

3.9540 

1185.0 

61.0 

167.2 

170.4 

OEC 

10 

73 

11 : 37 : 3 c 

323.7 

-21.5 - 

3.9944 — 

4105.0' — 



167.3 — - 

1-70*4 

DEC 

10 

73 

i 1 :ai : 55 

323.7 

-21.5 

3.9947 

1185.0 

61.0 

167.4 

178.4 

OEC 

1 0 

73 

11 :46i2C 

323.7 

- — 21.5 

3.9-549 - 

1 10S.fr 

61 .O 

467*6 - - 

170.fr- 

OEC 

1C 

73 

11 : 50:45 

323.7 

-21.5 

3.9550 

1185.0 

61.0 

167.6 

178.4 

OEC 

10 

73 

11 : 55 : 1 c 

323.7 

-21.5 

3.9951 

1 IBS * 0 

-64 ■ O 

467.6 - 

470,4 

OEC 

10 

73 

11 : 59:35 

323.7 

-21.5 

3.9951 

1185.0 

61.0 

167.7 

178.4 

0E£- 

-1 0 

7 3 

12 * 4 : c 

, 3-23V7- 

.... . ^21.5 

3 ^ g Q 5 . 1 — - * 

. , ... . 


. ... 1 Af “f 

4 TO A 

DEC 

1 0 

73 

12 : 8:25 

323.7 

-21.5 

3.9550 

1 185.0 

60.9 

— 

167.7 

t~r g m 4 - 

178.4 

DEC 

10 

73 

12 : 12:50 

323.7 

ai. fr - - - - 

- 3*9 9*9 

1 185.9 

-6fr*9 

467V? 

-470*4 

DEC 

10 

73 

1 2 : 1 7 : i £ 

223.7 

- 21.5 

3.9947 

1185.0 

60.9 

167.7 

178.4 

DEC 

10 

73 

12 : 21 : 4 c 

3 23.7 

- — 21. 5 

-3*994fr- 

-1 1 85. 0 

60.9 

167V6 

178.4 

DEC 

10 

73 

12 : 26 : 5 

223.7 

-21 .5 

3.9940 

1185.0 

60.9 

167.6 

178.4 

OEC 

1 0 

73 

12 : 30:30 — 

3E-3-. 7 - 

— — 31^5 

3 fjr . _ . 

Q « « „ 

6 0 9 

1 A f c 

4 *t0 A 

OEC 

10 

73 

12:34:55 

223.7 

-21.5 

3.9930 

1105.0 

60.9 

167«4 

178.4 

OEC 

1 0 

73 

12 : 39:20 

3 23.7 

-21.fr 

3*99-22 ■ - 

4-185.0 - 

- 60.9 

— 167.3 

— 170-*4 - 

DEC 

10 

73 

12:43:45 

223.7 

-21 .5 

3.9914 

1 1 65.0 

60.9 

167.3 

178.4 

OEC 

10 

73 

12 : 40:10 

- - 323*7 

-31*5, 

frfr00 3 

1 1-85.0- 

- 60.9 

Ifr7.fr 

- fr-7fr*4 

DEC 

to 

73 

12:52:35 

223.7 

-21.5 ‘ 

3.9692 

1105.0 

60.9 

167. S 

178.4 


1 2 1 67 0 

3 £3+7 

— 21 *5 

3 9 081 

L 1 66 .0 

- 60 a 9 

a t 0 

1 TO A 

DEC 

10 

73 

12 : 1:25 

323.7 

-21 .5 

3.5673 

1185. 0 

60.9 

168.1 

178.4 

-OEC 

1-0- -73 

13: s:so 

^3 34^0 

— -21.5 

- 3.9670 - 

1164.6 

- -6U1 

167.4 

~ 4J0*4 

OEC 

10 

73 

12 : 10:15 

324.2 

-21.4 

3.9074 

1101.4 

61*3 

167.5 

177.9 



B-8 


-C A VrifV«©3— ^49*1*7? Pne©«?©n—S9rt®V“«tt2“- -*“"&&<=— Sr& - - “ -VS“ 


PH6Wtf— 


0AT2 

T 5C5£ 

ASCENSION 

DECL2SSATEC* 

sp:^ rate 

wcwewTue: 

GAM?4A 

lakoa 

pmra&H 







r> ajJ f- t\ 

TT„M, J S^jT 



V* _* 

&&C 

-frO 

V-2- 

51 3 : n ^ S 4 0 

- -32 A~c>-£ — - 

— — - 




- --- -6-u — 




0£C 

£0 

7 3 

£3:i9: 5 

3 24 0 2 

-21 .4 

3.9693 

1 tsi.4 

Slo3 

E67 w2 

177*9 

-SEC 

■ s o 

73 

£3-323130 — 

-• 32^2 — 

- 22ro4- — ■ 


*- — 4 — — 

(S-Sr-vS* — 



- 

DEC 

£0 

73 

£3127:55 

324.2 

-21 .4 

3.9915 

1 161 .4 

61.3 

167 *t 

177*9 











DEC tO 

73 

13:36:45 

324„2 

-’21.4 

3.9931 

1191.4 

61.3 

167.1 

177.9 

DEC 

10 

73 

13:41 : 10 

-324.2 ~ 

24-*-4 

3^9936 

. * . | ^ 

- 64*3- 

“ 167. 2- 

477.9 

DEC 

to 

73 

13:45:35 

324.2 

-21.4 

3.9941 

1181.4 

61.3 

167.3 

177.9 

DEC 

10 

73 

135 50: c 

3 24.2 

-21.4 

-- 3.9944 

1 16 1*4 

61*3 

- -467 * 4- 

- 477 *9 

DEC 

10 

73 

13254:25 

324.2 

-21.4 

3.9947 

1181.4 

61.3 

167.5 

1 77.9 

-DEC 

10--73 

1 3258 150 

- 324.2- 



3.9940 

1 181.* 

6**8 

467*6-— - 

- 477.5- 

DEC 

10 

73 

14: 3 : 15 

324.2 

-21.4 

3.9950 

1181.4 

61*3 

167.6 

177.9 

DEC 

10 

73 

14: 7:40 

324.2 

-21-.* 

- -3. 9954 

110 4.4 

6-1.3- 

167.6 

177*9- 

DEC 

1 c 

73 

14:12: 5 

324.2 

-21.4 

3.9951 

1181.4 

61.3 

167.7 

177.9 

DEC 

10 

73 

145 165 3C 

324.2 

-31.4 

3.9951 

1 181.4 

61*3 

167.7 

477.9 

DEC 

1C 

73 

14:20:55 

324 .2 

-21.4 

3.9950 

1161.4 

61.3 

167.7 

177.9 

DEC 

1 0 

73 

14:25:20 

324.2 

- 2 1 *5 

3.9949 — 

1181.4 

61*3 

’ 467*7 - 

477.9 

DEC 

1 0 

73 

14:29:45 

224.2 

-21.4 

3*9947 

1 181.4 

61.3 

167*7 

177.9 

DEC 

10 

73 

14:34:1c 

324.2 - 

-21-4 

3.9-9 44- 

1 181 .4 

61.3 

167.6 

477.9 

DEC 

1 0 

73 

14238: 35 

324.2 

-2 1.4 

3.9940 

1181.4 

61*3 

167.6 

177.9 

DEC 

10 

73 

14143: C 

324.2 

-21.4 

-3*9935 

1181.4 

*4*3 

- 167.5 

177.9 

DEC 

10 

73 

I 424722S 

224.2 

-21.4 

3.9929 

1181.4 

61.3 

167.5 

177.9 

DEC 

10 

73 

14251 : 50 

324.2 

-31.4 

3.99 2 2-- 

1181*4- 

—6 1.3 - 

45A*4 

177,9 

DEC 

10 

73 

14256:15 

224.2 

-21.4 

3.9913 

1181.4 

6 1.3 

167.4 

177.8 

DEC 

ia 

73 

152 0240 

— 324.2 

— 31.4- 

— -3-*9 9 02 

1181.4 

64.3 

— 167-.4 

477.8 

DEC 

1 c 

73 

IS: 5: 5 

224.2 

— 21.4 

3.9691 

1 181 . 4 

61.3 

167.6 

177.8 

OEC 

10 

73 

15: 9:3C 

324.2 

-21-.* — 3. <5 680 

4^181.4 

6+-*3 - 

167.8- - 

- -- 1T7.8— 

DEC 

1C 

73 

15213:55 

324.2 

-21.4 

3.9673 

1 181.4 

61.3 

168.0 

177.8 

DEC 

10 

73 

15:18:20 

324.5 

'—24 « 3 — —3.0870 -- 

*479.7 

—6*. 6 

467^4 

477.6— 

DEC 

ia 

73 

15:22:45 

324.6 

-21.3 

3 • 9675 

1178.0 

61.7 

167.5 

177.4 

DEC 

10 

73 

15227: 10 

324.6 

-2 1 .3 — 

- 3.9564 

-4436.0- ... 

— 6 4*7- 

-- 16-7-.4 

477.4 

DEC 

10 

73 

15231 : 35 

324.6 

-21.3 

3.9695 

1178.0 

61.7 

167.3 

177.4 

OEC 

10 

73 

15236; C 

324.6 

-21.3 

3*9906 

1178.0 

61.7 

167.2 

477*4 

DEC 

1 0 

73 

15240:25 

3 24.6 

-21.3 

3.9916 

1178.0 

61.7 

167.2 

177.3 

OEC 

l a 

7 3 

15:44 : 50 

2 2-4. 6 

’“21*3 

3.9924 

H7A.rt 

6 1 4. 7 

. r -167 „ 2 ... , 

XT7*3 

DEC 

10 

73 

15:49: 15 

324.6 

-21.3 

3.9931 

1178.0 

61.7 

167.2 

177.3 

OEC 

10 

73 

15:S3:4C 

324.6 

-21*3- 

. 3*9937 

4178.0 

64.7- 

167.3 

477*3 

OEC 

1 0 

73 

is:£8: 5 

324.6 

-21.3 

3.9941 

1176.0 

61 .7 

167.4 

177.3 

OEC 

10 

73 

is: 2:3c 

324.6 

-21.3 - - 

,„*.J594S- 

- - 4 1 78 . 0 — 

64*6- 

. 4671*4-- 

477*3— 

DEC 

10 

73 

162 6255 

3 24.6 

-21.3 

3.9947 

1178.0 

61 .6 

167.5 

177.3 

dec 

io 

73 

16211:20 

324.6 - 

-74 .3 

3. 0949 

1176*4) 

64*6- 

457*6 — 

177*3- 

OEC 

10 

73 

16215:45 

3 24.6 

-21.3 

3.9950 

l 178.0 

61.6 

167.6 

177.3 

OEC 

10 

73 

16:20:i0 

3 2-4.6 

-21^*3— 

3*9954 

*1-78* 0 t - 

54.6- 

467.6- 

477*3- 

OEC 

10 

73 

16224:35 

324.6 

-21.3 

3.9951 

1178.0 

61.6 

167. T 

177.3 

-OEC 

ID 

73 

16229: 0 

324.6 

-24.3 - 

3*9954 

1 178*0- 

-61.6 

167*7 - - 

137*3 

OEC 

10 

73 

16:33725 

324.6 

-21.3 

3.9950 

1178*0 

61.6 

167.7 

177.3 

l\Ci* 1/1 ^ <l 

1 A * 

334,6 

—31,3 

1 ■ 9 94 R 

1 t 7Q . 0 

— - - 6 1.6 

167.7 

1 77 t 5 

DEC 

10 

73 

w 

1 ** 

♦♦ 

<r 

2 

224.6 

-21.3 

3. 9 946 

1178*0 

61.6 

167.7 

177.3 

1 <4 

73 

16:4*: 4« 

3 24.6 

- . -21 *3 

__3*9*43 

1178.0- - - 

51*6 

167.6 

177*3 

DEC 

10 

73 

162512 5 

324.6 

-21.3 

3.9939 

1178.0 

61.6 

167.6 

177.3 



**-■ «• - c — ATTtTtru* - ttnwmmr 


PPgPICTG B S T A RT TI ME 


OATE 


RKttT 

ASCENSION 

< OE< 5R EC5» 


DECLINATION 


SPIN RATE 
H*p» * 3 -»- 


MCHENTUM 


GAMMA 

<DC6il 


LAND A 
< PEG . 


PH I MAX 


w 

I 

CD 


- DEC 16- T 3 — 
DEC 10 73 
^ 1-0—73 - 
DEC 10 73 
-DEC 10^73 
OEC 10 73 
DEC lO 73 
OEC 10 73 
DEC 10 73 
OEC 10 73 
OEC 10 73 
DEC 10 73 
DEC 1073 
DEC 10 73 
DEC 10 73 
DEC 10 73 
“DEC— 1 O 73 
OEC 10 73 
DEC 10 73 
DEC 10 73 
DEC 10 73 
DEC 1C 73 
-DEG- tO 73 


16 “55 1 30 — 

16:59 :ss 

17+-*: ED 
17: e:45 
17+I3 H-* - 
17:i7!35 
17122! O 
171 26 : 25 

it: ao isc 
17:35: 15 
17:39140 
17:44: 5 
17:49! 3C 
17152155 
17:57:20 
161 1145 

1 e: 10 

19110135 

ie:is: 0 
10:19:25 
10:23:50 
10:2s: 15 
10:22:4c 


- 204*6 
224*6 
-2 24i+r 


-21*3 

2W-3- 


3**966 
3*9926 
-3*99*1— - 


1 - 1 - 7 3 >-6- 
1 170,0 
-1173*0 


61*6 

61*6 

-6* * 6 


t6 T« S- 
167*5 
1 6 7 * 5 


1 7 7*- 3 

177*3 

+77i-3- 


DEC 10 73 
DEC 10 73 
DEC 10 73 
DEC 10 73 


19: 9: 0 
19:12:25 

19:16:50 

I9:ai : is 


DEC- 1-6 73 — 15125:40 


DEC 10 73 


is:3o: s 


OEC 

10 

73 

19 

:34; 30 

DEC 

10 

73 

IS 

: 38:55 

OEC 

10 

73 

IS 

:4-3 : 20 

OEC 

1 0 

73 

19 

:47 :45 

DEC 

10 

73 

154 52 MO 

DEC 

10 

73 

IS 

: 56 : 35 

DEC 

1 0 

73 

20 

1 l x 0 

DEC 

10 

73 

20 

: 5:25 

DEC 

10 

73 

20 

: 9:50 

DEC 

10 

73 

20 

M4J 15 


324*6 



-21,3 

3*9911 

1170*0 

61*6 

167*4 

177*3 






167 ttS 

ITTirf- 

224*6 

-21*3 

3,9009 

1170.0 

61*6 

167.6 

177*3 

3E4--6-- - 

S+v3 

3^90 79 -- 

— 1 -1-76 , 0- — 

6-1 .6“ 

167.8- 

— 1 77 . 3 

224-6 

-21,3 

3*9072 

1 170.0 

61.6 

166*0 

177.2 

225*0 

— 21*2 

3^9 e7-i 

n74.e 

62.0 — 

167*4 

176*9 

325,1 

-21*2 

3-9675 

1 174*0 

62.0 

167*5 

176*6 

3 25-* 1 

— -gj ig 






325,1 

-21*2 

3*9696 

1 tmiu 
1174.0 

Or* w — 
62.0 

167*4 

176.6 

-5 25*1 

—21 *2 - 

5*590*- - 

1-+74*G - 

6 2-. 0 

167 tt3-— - 

1 76* 8 

225,1 

-21,2 

3. 9917 

1 174.0 

62.0 

167.3 

176.6 

225*1- - 

-«b2 

3.5525 

1174*0 

62 rO™ ■- — 

1 67- *3 

+76*9 

225,1 

-21.2 

3*9932 

1174.0 

6 2-0 

167*3 

176*8 



DEC 10 

73 

16:37: 5 

3 25.1 

-21.2 

3*9951 

DEC 10 

73 

10:41 : 30 

325.1 

-2 1 ^2* — 

— -3-995+ 

DEC 10 

73 

10:45:55 

325.1 

-21.2 
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APPENDIX C - OPTIMAL ATTITUDE CONTROL 


COMPUTER SIMULATION FOR TIME OPTIMAL OR ENERGY OPTIMAL ATTITUDE CONTROL ■ 
OF SPIN-STABILIZED SPACECRAFT, 


Robert D. Woolley 

(Computer Sciences Corporation, silver Spring, Md.) 
Roger D. Working 

(NASA-Goddard Space Flight Center, Greenbelt, Kd.) 


SUMMARY 

An original technique for determining the optimal magnetic 
torque strategy for control of the attitude of spin stabilized 
spacecraft is presented. By employing Lagrange multipliers and 
the Calculus of Variations, optimal control equations are de- 
rived which define minimum time and minimum energy attitude 
maneuvers. Computer program algorithms to numerically solve 
these optimal control equations are also described. The per- 
formance of this technique is compared with a commonly employed 
planning method. 

INTRODUCTION 


For several years, the orientation of most near-earth 
spacecraft has been controlled by the interactions of magnetic 
fields . Figure 1 illustrates a controllable electromagnetic 
coil whose axis is parallel to the axis of rotation (spin axis) 
of the spacecraft on which it is mounted. When electric cur- 
rent flows through this coil, the resultant magnetism interacts 
with the earth’s magnetic field to produce a torque on the 
Spacecraft, By controlling the magnitude and polarity of the 
electric current, controlled torques are produced which change 
the spin axis orientation. 



SPACECRAFT'S 



Figure 2. QOMAC Coil Switching Technique 

This paper presents a technique to mechanize the control 
planning process using high-speed computers and optimal con- 
trol theory. The computer program developed for this tech- 
nique will be used to support the Atmosphere Explorer (AE) and 
Orbiting Solar observatory (OSO) a a tel li tea. It will minimize 
either the coil energy or the elapsed time required to complete 
maneuvers between specified attitudes. The option selected de- 
pends on the immediate mission requirements, 

THE ATTITUDE STATE EQUATION 

Consider a spin stabilized spacecraft system containing 
a nutation damping system. This system maintains collinearity 
between the spacecraft’s angular velocity vector and one of its 
principal body axes. Defining the spin axis as their ccranon 
direction, the angular momentum is 


f-um ^ i * h ° raa 9 ni tvde of the spacecraft’s angular momen- 
tum, and s is a unit vector in the direction or the spin axis, 
expressed in an inertial coordinate system. 


Figure 1, A Spin-Stabilized Spacecraft Employing 
Magnetic Attitude Control 

In early missions, commanded changes of this electric cur- 
rent were limited by the necessity for ground station contacts 
at the times of change. The magnetic control systems were left 
in a given state for a period of orbits or daya. During this 
time, the spin axis would precess at a very alow average rate. 
These early control systems were intended to maintain spacecraft 
sun angles and compensate for unavoidable environmental dis- 
turbance torques. 


t „. ^Pin-stabilized spacecraft utilising magnetic .ttitude con- 

"V ? “ d built with tha 4)118 ° f «* onboard bidiree- 
tlonal electromagnet aligned to the spacecraft spin axia. The 
?, r J. y 0f , te resulting magnetic moment la determined by the 
.f' C ‘ 0 " of electrical current through the coil, the magnitude 

rent 18 deter) " ined by the magnitude of coil cur- 

* . * Theae design characteristics are summarized by the fol- 
lowing equations. ’ 


( 2 ) 


The launch of TIROS IX introduced an innovation in magnetic 
control system design, it was known that a more efficient con- 
trol system would switch the coil polarity four times per 
orbit (1J, [2J, 13], This system, depicted in Figure 2, was 
called Quarter Orbit Magnetic Attitude Control {QOMAC) . The 
quarter orbit switching was accomplished by an onboard timer 
which reduced the reliance on ground station contact. The QOMAC 
manouver continued as long as was necessary to achieve the de- 
sired spin axis precession. Control capabilities increased 
dramatically from a few degrees per week to about five degrees 
per orbit. 

Recently designed spacecraft are being equipped with com- 
mand memories. With this capability, the exact time at which a 
coil state should change can be preloaded into the command 
memory, This has added tremendous flexibility to the use of 
magnetic control systems and has improved the efficiency of 
these systems even further. This is especially true for 
eccentric orbits and low inclination orbits where QOMAC is not 
readily adaptable, with this increased flexibility, control 
planning processes become more complicated. 


where M ie the electromagnet • s magnetic maoent, m is the 

of “ 49n ® tic <po»eibly a function 

state. * u i. tha normalrted remotely commanded coil 


, ZZ ’ w 8 P«=«ciart resulting from the 
n .Sc fis“ L d * IeCtrcm “ 9n8 t «tth the earth's mag- 


-magnetic — _ _ 




rn.9n.tir flit deni?t t ►r Pr ? Mnt<,tlC "* ° f the 1 °“ 1 

n«ar!£rtf “ ^ ° the / lclnlty 01 the spacecraft. In 

thesis ’ ma9netic fieid 1« approximately equal to 

“li fi * ld - Thi8 18 m. themati oally character- 

in1 a « 8 *:f f l'^r t / CCUr#Cy bV * >Pherlsal h armoni ^expansion 

' “ d altitude HI. [5] ? Thu., 

magnetic flux densitu u nea *-*arth orbit, the local external 
tiS alon“ y ' ~ b * obt » ln8d •• 8 function of 
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The spin ax is state equation is governed by Newton's 
second law applied to the rotational motion of the spacecraft 
in inertial space. Thus 

*= £ * mBu + T^ (5) 

where is the vector sum of environmental disturbance 

torques acting on the spacecraft. Significant disturh.inr'* 
torques are caused by the interaction of the residual magnetic 
moment of the spacecraft with the earth's magnetic field, and 
by offsets between the spacecraft center-of-mass and the lines 
of action of solar pressure, aerodynamic drag pressure, and 
gravity. These offsets and residual magnetic moments tend to 
rotate with the spacecraft, so that their effective average 
direction is approximately aligned with the spin axis. The 
disturbance torques can thus be approximated by 

Id " 1 x H (6) 

where JD in the sum of the average generalized forces multi- 
plied by their offsets along the spin axis. 

Spin stabilized spacecraft generally have extremely small 
variations in the magnitude of the angular momentum vector. 

Ihia is partial iy because spin axis components of environmental 
torques tend to cancel themselves effectively during each spin 
rotation po>:iod, and partially because the angular momentum 
magnitude is generally quite large. Since angular momentum mag- 
nitude is maintained within a close tolerance of a nominal value 
by periodic actions of a different onboard torquing ay stem, it 
is considered to be a constant or a known time function. Thus 
wa have tho following spin axi3 attitude state equation! 



For most spin-stabilised spacecraft, the coil controlling on- 
board electronics do not allow continuous variation of coil 
current magnitude. Instead, the coil current nay be off or 
on with either polarity at a designed value. For these space- 
craft, Equation (3> is roplaced by 

u ** -1, 0, or <*1 (8) 

THE CONTROL PLANNING PROBLEM 


THE OPTIMAL CONTROL EQUATIONS 

in view of Equation {8} , the amount of energy used by the 
coil for a given maneuver is proportional to the functional 


r 


iui dt 


( 9 ) 


The elapsed time for a maneuver is the functional 




Using the definitions 


< 10 > 


^ " iui (11) 

l 2 ° 1 (12) 

the performance functions which arc to be separately minimized 
may be compactly expressed by the single integral equation 


-.A 


dt 


( 13 ; 


The quantity being minimized, J , is augmented by the state 
and control equations as follows j 


J i ■ j [ L i + i T ^. * [| £ u * ff S.] - dt)* “<“ 2 - !>] 


dt (Id) 


where X and u are Lagrangu multipliers with 


[>0i£u -1^0 


0 if u - 1 < 0 


( 15 ) 


The attitude state equation, Equation (7), is time varying, 
nonlinear, and not totally controllable. The time-varying 
nature of the magnetic and disturbance vectors, £ and D , pre- 
clude the use of classical controls techniques which employ fre- 
quency dcoain analyoio. Nonlinear effects arise because the 
control variable, u , and the attitude, s , are multiplied to- 
gether in Equation (7). The attitude is not totally controllable 
because instantaneous disturbances and the desired attitude 
motion possess two degrees of freedom while the instantaneous 
coemondable torque has only one. Therefore, it is theoretically 
impossible to precisely follow many continuously specified atti- 
tude trajectories , The addition of the control constraints in 
either Equation (3) or Equation (8) to Equation (7) almost 
guarantees the fruitleasneas of an intuitive approach to con- 
trol planning. 

Attitude control pliinning may be accomplished with time 
doraain analysis if seme additional constraint is placed on the 
control beyond those of Equationo (7) and (8J . For instance, 
tho QCMAC method constrains the coil to remain energized 
throughout an active control period, alternating polarity at 
proctso quarter-orbit intervals. Thus, the QOMAC control prob- 
lem is reduced to determining the desired duration of the active 
coatrol period and the constant relative phasing between the 
orbit and coil alternations. 

An optimal control lav adds the constraint that a physi- 
cally significant cost function is minimized. The objective of 
the control technique to be presented is to determine the dis- 
crete utt) control history which causes the spin axis atti- 
tude to change from to £2 in an optimal manner. The 

direction, £<t) , must alwayB Batisfy Equation (7), and the 
sense of optimality is either minimum time or minimum electri- 
cal energy to activate the electramagnat. Minimum time 
maneuvers achieve greater attitude motion during a limited 
period than any other maneuvers. Hence, rapid attitude changes 
are acctxnplished moot effectively. Minimum energy maneuvers 
achieve a desired final attitude in a specified time period, 
with as little coil energy expenditure as possible. 


Equation (14) is integrated by parts and its first variation 
eat to zero as a necessary condition for minimization of J. . 
Since the endpoint conditions are fixed, and the Lagrange i 
multipliers are selected ao that the remaining terns in the 
first variation of with respect to variations in the con- 

trol, u , and the attitude, , are zero, any optimal maneuver 
minimizing Equation (14) must satisfy the following Euler- 
Lagrange vector equations [6Js 


f?(h * ± T [i * (ff* * kzb (St i*)i* - »j ■ 0 US, 

L. j tl + * T [£ « (f £u ♦ i d)]+ (fj x T )i ♦ U |u J - 1)| = o t U7| 


The Lagrange multipliers, the attitude, and the disturbance 
torques ara not explicitly functions of the coil state, u 1 
therefore, partial derivatives of ^ , \i , £ , and d with re- 
spect to u are zero. Coil control histories for minimum 
energy and minimum time maneuvers are defined as u. (t) and 
u^<t) , respectively. Equations (11), (12), (15), and (16) arc 

ccjfibined to produce 


L ’ 

<s * B)j 

(IB) 

. - - S «[|i T 

(s x B)J 

(19) 

the decdzone 

function 


/ +1 if 1 

< X 


;) - < 0 if -1 

£ x £ +1 

(20) 

V-l if x 

< -1 
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and where SGN represents the signum function 

/+1 if o < x 

SGN(x) ** < 0 if x « 0 (21) 

w if x < 0 

The variables , 1^ , u , and U are not explicitly functions 
of the attitude; therefore , their partial derivatives with re- 
spect to £ are zero. For either the minimum time or the mini- 
m-urn energy optimal maneuvers. Equation (17) is reduced to 



Equations (7), (18), and (22) or Equations (7), (19), and (22) 
define minimum energy or minimum time optimal maneuvers, respec- 
tively, in the presence of any normal disturbance torques. How- 
ever, for certain types of disturbance torques, such as those 
caused by residual dipole moment or aerodynamic forces, the 
vector £ is approximately a function of time only, independent 
of ♦ For these types of disturbance torques. Equation (22) 
becomes 


COMPUTER ALGORITHMS 

The optimal control equations have two undetermined 
parameters which must be correctly selected in order to 
achieve the final attitude specified by Equation (30) , One 
of the parameters appears explicitly in Equation (31) , the 
other parameter is implicit in Equation (28) . TO express this 
other parameter in equation form, a right-handed system is de- 
fined about the initial attitude. Introducing two mutually 
perpendicular unit vectors, and , which are perpendic- 
ular to the initial attitude vector, i.e., 

" 0 

S(t x ) “ | 2 ) - +1 (33) 

Equation (28) can be replaced by 

2/t^) - a COS 4> + a 2 sin $ , o s 4 < 2* (34) 

The two undetermined parameters are then C and 4 for the 
minimum energy case, and t f and 4 for the minimum elapsed 
time case. 


iti’ i* [|i“ ♦ iU] < 23 > 

Equations (7) and (23) maintain both a constant magnitude of 
1. and a constant angle between X_ and £ . Since the com- 
ponent of ^ along £ does not affect u through either 
Equation (18) or (19), that component may be arbitrarily set to 
zero. The resulting X_ has a magnitude and a direction per- 
pendicular to £ which define 

c - — == ( 24 ) 

JK 

A ■ £ x i c < 2S > 

By the definition of the signum function given in Equation (2i> , 
the coil control for a minimtsn time optimal maneuver is not in- 
fluenced by the magnitude of , only by the direction of . 
Furthermore, as the magnitude of approaches infinity, Equa- 
tions (18) and (19) become equivalent. With the above consid- 
erations and a specification of required initial and final 
attitudes. Equations (24) and (25) are substituted into Equa- 
tions (18) and (23). The following resultant optimal equations 
define both minimum time and minimum energy optimum maneuvers! 


d a a I'm 11 

[iEu.-Dj (7, 

it i- i* [flu + !d] ( 26 ) 

1 +1 if C < | B*£ and t s t f \ 

0 if -C * £ scor t { < t > (27) 

-1 if * -C and t s t f / 

" 0 * 2e) 

t>l , £(t^) - specified (29) 


Because the boundary conditions of Equations (29) and (30) 
are split between the initial and final endpoints, no explicit 
equations are employed to determine the parameters 4 and c 
or t f , Instead, Equations (7) , (26) through (29) , and (31) 
are successively solved by a fast simulation routine. A two- 
dimensional numerical search algorithm drives this routine, 
varying values of the appropriate parameters on successive 
simulations to converge on the final attitude specified by 
Equation (30). Computer time economy dictates that both the 
fast simulation and the two-dimensional search be as efficient 
as possible. 


The commendable magnetic precession vector appearing in 
Equation (27) is a function of time alone, independent of 
attitude. The geomagnetic flux density vector, £ , is a known 
function of the spacecraft's position with respect to the rotat- 
ing earth; the spacecraft's inertial posit ion- time dependence 
is given by non-Keplerian formulas which include effects of 
terrestial oblateness and atmospheric drag. The commendable 
magnetic moment, m , is a function of the voltage available to 
power the onboard electromagnet. This voltage may have a dif- 
ferent value when the spacecraft's position is in the earth's 
shadow than when the spacecraft's position is in sunlight, 
because of differing characteristics of batteries and solar 
arrays. The comroandable magnetic moment may also be set to 
zero during any portion of the period of interest in order to 
mechanize requirements that the coil be off during those 
periods (see Equation 27) . 


Although the command able magnetic precession vectors 
could be recomputed during each execution of the fast simula- 
tion routine, much computer time is saved (at the expense of 
some storage space) by initially filling a 3-by-N array with 
N computed 3-by-l commandable precession vectors for each N 
integration time step spanning the period of interest. The 
fast simulation routine ig further streamlined by the use of a 
direct precession vector integration scheme. Equations (7) and 
(26) may be rewritten as 


at 1 ’ 3 * [ir ! u * h s] + ( i * 2> * [f i u + £ e] os) 

3t a • % * [| 2“ ♦ £ e] + <i - a> » [| £u + ! d] (3«) 


** * i‘V ’ k 


specified 


(30) 


where £ and ^ may be arbitrary. For a constant coil con- 
trol period where 


tj ■ t^ , C i 0 if minimum energy 

(31) 

C ■ 0 » s t { 1 t 2 if minimum time 

Since Equations (7) and (26) have identical form, they can both 
be numerically solved in one attitude simulation program. 
Equation (29) fixes the initial attitude s(t ) and Equa- 
tion (28) constrains the initial pseudo-attitude £{t ) to be 
perpendicular to £(t.) • Equation (27) gives the algorithm 
-used for computing the optimal coll control throughout the simu- 
lation . 


they ore selected 


U J ' * t * Vi 

(37) 

1 * itty 

(38) 


(39) 
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Equation® (35) through (39) are combined, neglecting the small 
terms in Equations (35) and (36) , and integrated to yield 



Analysis of the approximations, Equations (40) and (41), shows 
that the directional approximation errors can be kept extremely 
email while maintaining appreciable time steps between t and 
t j . For small motions through an angle (radians), the 

per unit directional error introduced by Equations (40) and (41) 
is 


(42, 

For example, a 90-degree attitude maneuver simulated using these 
approximations and Aifi as large as one degree would be in error 
by about one hundredth of a degree. Of course, each applica- 
tion of Equations (40) and (41) is followed by unitizing the 
results to form £(t^ +1 J , ${t^ +1 ) . 

Since the integrals in Equations (40) and (41J do not in- 
volve the attitude, it is not necessary that they be recomputed 
in each execution of the fast simulation program. Much computer 
time ig saved by initially filling a pair of 3-by-(N+l) arrays 
with the (N+l) computed 3 -by-l integral contaandable precession 
vectors and ths (N+l) computed 3-by-l integral disturbance pre- 
cession vectors for each of the N+l integration time points 
spanning the period of interest. These integral commendable and 
disturbance precession vectors are directly computed as the 
Riemann integral of the commendable and disturbance precession 
vectors. Position and time dependent parameters essential to 
compute the canmandable and disturbance precession vectors are 
provided by a call to an external routine. 

The fast simulation program steps through every entry of 
its cocuandable precession vector table, using Equation (27) 
in every case to determine the instantaneous commanded coil 
state. Equations (40) and (41) are applied every M entries 
to avoid excessive approximation error. Whenever a coil state 
change is detected, a search is instituted using Equations (41) 
and (27) to determine the precise switching time. After it is 
located, both Equations (40) and (41) ace applied before con- 
tinuing to step through the table, checking Equation (27), 

The two-dimensional numerical search algorithm which calls 
the fast simulation program is baaed on a variable separation 
technique. A right-handed system is defined about the final 
attitude resulting from no control, such that 

S(t > « <£, « 6 2 > - +i («> 

CONTROL 

£ x x fi 2 - 0 (44) 

then for any other final attitude, s^Uj), can define angles 
P k and & k such that 

Lit J - £<t ) cos p ♦ ^ cos G k sin 

d NO CONTROL 

+ fij ain 0 k sin p k 


For both- the minimum elapsed time case and the minimm coil time 
case, it is assumed that 8 is monotonically related to 4 
We also assume that p is monotonically related to C in the 
minimum energy case, and to tf In the minimum time case. The 
search technique, making use of these assumptions tries to ad- 
just * and either C or t f in order to drive 6 and p 
to the desired 6 and p . The adjustments are made in em- 
bedded bracket-method, one-dimensional searches. These 
tions have been justified by the convergence of the search 
algorithm for all tested cases. 

The inner search algorithm adjusts only the value of C 
(or t. ) , utilising in its calls to the fast simulation rou- 
tine the $ passed to it from the outer search. It initially 
tries C - 0 (or t f ■ t 2 ) to establish whether or not an 
angular motion of PnzsUEV be achieved. If the resultant 
p confuted by the fast simulation routine is less than the 
desired p , control is passed back to the outer search, if, 
however, the desired magnitude of motion can be achieved, the 
inner search successively increases C (or decreases t f ) 
until a p less than the desired p results. Having bracketed 
the desired amount of motion with two values of C (or tf ) , 
the inner search iteratively converges on a value producing 
the desired motion by repetitively testing values in the inter- 
val and letting them replace the appropropriate upper or lower 
bracket value to reduce the bracket size. Control is passed 
back to the outer search algorithm after the bracket size is 
reduced to the computer's numerical precision size. 

The outer search algorithm starts with an initial esti- 
mate of the correct value of $ , then passes control to the 
inner search algorithm. After the inner search algorithm has 
returned control to the outer search, 0 is computed and com- 
pared with the desired 0 • Then 4 is either increased or 
decreased in repetitive calls to the inner search, until the 
outer search has established a pair of *)-yalues with associated 
0- values bracketing the desired 9 . After bracketing is 
achieved, convergence proceeds similar to the inner search, by 
repetitively subdividing the interval, down to the numerical 
precision of the computer. Upon convergence of the outer 
search, the desired attitude maneuver has been obtained. The 
coil state history is that determined during the last execu- 
tion of the fast simulation routine. 

The average number of fast simulations executed to achievp 
ultimate convergence is approximately the product of the 
average number of bracket subdivisions performed in each of the 
two one- dimensional searches. Thus, the required computer time 
is quite sensitive to the method of subdividing the intervals. 
The conservative approach of interval halving requires 25 steps 
per one- dimensional search. Use of this method requires about 
625 fast stmulaUona for ultimate convergence. A linear in- 
terpolation approach can complete a one-dimensional search in 
as few as 4 steps if the linearity assumption is approximately 
valid, requiring about 16 fast simulations for ultimate con- 
vergence. Drastic departures from the assumed linearity can 
make this method even less efficient than the interval halving 
method. 

The bracket subdivision algorithm used in the developed 
program is a weighted average between linear interpolation and 
Interval halving. Weighting factors are redetermined prior to 
each iteration on the basis of the departure from linearity de- 
tected in the previous iteration. This self-adaptive feature 
results in consistent achievement of the ultimate convergence 
within 20 to 50 fast simulations. 

If the desired final attitude is unachievable during the 
period of interest, this search scheme converges to a final 
attitude on a great circle connecting the desired final atti- 
tude with the M no control" final attitude, as close as possible 
to the desired final attitude. However, this program can also 
identify the set of achievable final attitudes when operating 
it in a different mode. 

The set of all final attitudes achievable during a control 
period is bounded by the set of all final attitude# achievable 
using maximal time optimal control throughout the period. 

Thus, by setting 


( 46 ) 

( 47 ) 
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and stepping $ from 0 to 360 degrees on successive executions 
of the fast simulation program, a contour passing through the 
successive attitudes encloses all possible final attitudes 
achievable during the control period utilizing any magnetic 
attitude control planning scheme. 

EVALUATION OF OPTIMAL MAGNETIC ATTITUDE CONTROL 

Practical application of optimal control to real attitude 
control planning tasks necessitates access to a computer. Al- 
though the time required for the computer solution is minimal, 
a manual numerical calculation of an optimal attitude maneuver 
would be prohibitively time consuming. This computer require- 
V ®snt compares unfavorably with previously employed planning 
methods. Approximate algebraic formulas exist for both QOMAC 
and continuous control maneuvers, allowing an analyst to 
manually plan maneuvers with few numerical calculations. But 
the benefits gained from optimal maneuvers warrant computer 
availability. 

Optimal magnetic attitude control has two main advan- 
tages, Coil utilization is more efficient for optimal control 
maneuvers than for any other type of maneuver. Also, the use 
of optimal control automates the planning process, thus re- 
ducing the analyst's involvement below the necessary levels of 
involvement characteristic of other planning methods [7], 

Coil utilization Is significantly more efficient for opti- 
mal control than for QOMAC, over a large class of orbits. The 
only cases where the improvements in efficiency are sometimes 
small are with low altitude polar orbits. In those cases, 
optimal maneuvers are similar or identical to QOMAC-gene rated 
maneuvers. 

Table 1 compares minimum time and minimum energy optimal 
maneuvers with a QOMAC maneuver. 

Table 1. Performance Comparison of Optimal and 
QOMAC Maneuvers 


MANEUVER TYPE 

MANEUVER 
COMPLETION TIME 

ENERGY 
(COIL ON-TIME I 

MINIMUM TIME OPTIMAL 

270 MINUTES 

270 MINUTES 

MINIMUM ENERGY OPTIMAL 

457 MINUTES 

174 MINUTES 

QOMAC 

467 MINUTES 

457 MINUTES 


The maneuvers are between identical attitudes and start 
at the same initial time in the same orbit (see Table 2) , 

Table 2. Parameter Values Used in Comparative Maneuvers 


PLANNING PARAMETERS 


INITIAL TIME 

DEC. 21. 1073, 0 HRS. 0 MIN. 0 SEC 

INITIAL RIGHT ASCENSION OF SPIN AMIS 

10 DEGREES 

INITIAL DECLINATION OF SPIN AXIS 

0 DEGREES 

FINAL RIGHT ASCENSION OF SPIN AXIS 

0 DEGREES 

FINAL DECLINATION OF SPIN AXIS 

0 DEGREES 

j SPACECRAFT PARAMETERS | 

comma* dable magnetic moment of coil 

08 AMPERE TURN -METER 2 

ANGULAR MOMENTUM ABOUT SPIN AXIS 

134.4 N EWTO N-ME TE R-SE C 

| ORBIT PARAMETERS j 

ECCENTRICITY 

0.002 

SEMIMAJOR AXIS 

0.700 KILOMETERS 

inclination 

OS DEGREES 

RIGHT ASCENSION OF ASCENDING NODE 

□ DEGREES 

ARGUMENT OF PERIGEE 

0 DEGREES 

MIAN ANOMALY 

□ degrees 

EPOCH 

DEC. 21. 1S73: 0 HRS, 0 MIN, 0SEC 


The chosen orbit is low altitude, nearly circular, and 
only 22 degrees away from a polar orbit, to present a case where 
QOMAC performs well. But the minimum time optimal maneuver 
(Figure 3) is still considerably faster than the QOMAC maneuver 
(Figure 5) . The minimum energy optimal maneuver (Figure 4) was 
computed by allowing the same amount of .time for maneuver com- 
pletion as required by QOMAC. But the minimwn energy optimal 
maneuver uses 60 percent less energy (coil-on- time) than the 
QOMAC maneuver. These improvements become considerably more 
dramatic with orbits of lower inclination, higher altitude, or 
greater eccentricity. 



Figure 3. Minimum Time Optimal Maneuver 



Figure 4. Minimum Energy Optimal Maneuver 
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Figure 5, QOMAC Maneuver 
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Numerous codifications to the QOMAC method have been de- 
veloped for special cases and used for control planning. These 
modifications attempt to improve the coil utilization efficiency 
vjith minimal change to the analytical approach of the basic 


QOMAC method. Of course, 
these pseudo-QOMAC methods 

P^OOn fnv .»k4 J-l* 

- wi u»ejf 


coil utilization efficiencies of 
, even when applied to the special 
developed I Con never surpass the 


coil efficiency of the optimal method. 


Optimal magnetic attitude control automates the planning 
process through its algorithms' form and flexibility. The op- 
timal control algorithms are orbit independent. They therefore 
apply directly to elliptical orbits. They can be used for a 
control period including only a portion of an orbit, or even 
for a control period which includes orbit adjustment maneuvers. 

The optimal- magnetic attitude control algorithms developed 
hare compute a maneuver from the input specified initial atti- 
tude to precisely the input specified final attitude, if that 
final attitude is within the region achievable during the 
planning period. Planning methods utilizing an integer num- 
ber of quarter-orbit coil pulses cannot, in general, precisely 
achieve a specified final attitude. The optimal control plan- 
ning method eliminates the need for a "fine tuning" terminal 
portion of a maneuver as required by a modified QOMAC tech- 
nique, thus reducing analyst requirements. 


The command able region boundary of achievable final atti- 
tudes may be automatically computed for any planning period, 
using the optimal control algorithms. An example of a com- 
mendable region boundary for an 6-hour control period using the 
parcaeters of Table 1 appears in Figure 6. None of the atti- 
tudau outside the boundary can be achieved during the control 
period using magnetic attitude control, no matter what method 
of planning is used. Any attitude inside the boundary can be 
achieved during the control period. Any attitude on the 
boundary itself can be achieved during the control period but 
only by employing optimal control. 



methods, between a control planning computer program, and an 
accurate attitude prediction (simulation) computer program. 
Elimination of the discrepancies reduces the need for analysis. 

The optimal magnetic attitude control algorithms have the 
additional flexibility that a control period may include 
several input-specified, coil-off periods. A control program 
operator, knowing of scheduled spacecraft data collection 
periods during which the coil must be kept off, simply speci- 
fies those periods to the program. The commandable region 
boundary, and any optimal maneuvers, are then automatically 
computed by the program subject to the constraint that the 
coil be kept off during the specified coil-off portions of 
a control period. 

By first presenting to an operator an accurate display of 
achievable final attitudes subject to any operator-specified 
time constraints, and then by generating a set of commands 
which precisely achieve the operator-selected final attitude, 
the optimal control algorithms effectively automate the con- 
trol planning process. The generated set of commands mini- 
mize coil usage over all possible maneuvers meeting the 
opera tor- specified conditions. 
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Figure 6. Commendable Region Boundary for an 
Eight-Hour Control Period 


By providing the capability of computing the commandable 
region boundary for any control period, the optimal control 
algorithms succinctly display the set of available options to 
a program operator. This flexibility greatly reduces the 
analytical expertise necessary to select a desired final atti- 
tude, or a control period, or both. 

Approximation errors in the optimal magnetic attitude con- 
trol algorithms can be kept arbitrarily small. Consequently, 
any degree of desired planning accuracy can be achieved. An 
accurate geomagnetic field model can be used, as opposed to 
the canted dipole model assumed by most other methods. Dif- 
ferent comaandable coil strengths may be used for the sunlit 
and dark portions of a planning period. The effects of any 
codolled disturbance torques can be directly included. These 
characteristics of the optimal control algorithms can elimi- 
noto the discrepancies, innate in other control planning 
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